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ABSTRACT 

We present cosmological parameter constraints based on the final nine-year WMAP data, in conjunc- 
tion with a number of additional cosmological data sets. The WMAP data alone, and in combination, 
continue to be remarkably well fit by a six-parameter ACDM model. When WMAP data are combined 
with measurements of the high-Z cosmic microwave background (CMB) anisotropy, the baryon acous- 
tic oscillation (BAO) scale, and the Hubble constant, the matter and energy densities, ^bh^ , ^ch^, 
and f^A, are each determined to a precision of ~1.5%. The amplitude of the primordial spectrum is 
measured to within 3%, and there is now evidence for a tilt in the primordial spectrum at the 5cr level, 
confirming the first detection of tilt based on the five-year WMAP data. At the end of the WMAP 
mission, the nine-year data decrease the allowable volume of the six-dimensional ACDM parameter 
space by a factor of 68,000 relative to pre- WMAP measurements. We investigate a number of data 
combinations and show that their ACDM parameter fits are consistent. New limits on deviations from 
the six-parameter model are presented, for example: the fractional contribution of tensor modes is 
limited to r < 0.13 (95% CL); the spatial curvature parameter is limited to fJfc = — 0.00271q;oo38; the 
summed mass of neutrinos is limited to '}2,m^ < 0.44 eV (95% CL); and the number of relativistic 
species is found to lie within A'off = 3.84 ±0.40, when the full data are analyzed. The joint constraint 
on TVcft and the primordial helium abundance, Ihc, agrees with the prediction of standard Big Bang 
nucleosynthesis. We compare recent Planck measurements of the Sunyaev-Zel'dovich effect with our 
seven-year measurements, and show their mutual agreement. Our analysis of the polarization pattern 
around temperature extrema is updated. This confirms a fundamental prediction of the standard 
cosmological model and provides a striking illustration of acoustic oscillations and adiabatic initial 
conditions in the early universe. 

Subject headings: cosmic microwave background, cosmology: observations, early universe, dark matter, 
space vehicles, space vehicles: instruments, instrumentation: detectors, telescopes 



1. INTRODUCTION 

Measurements of temperature and polarization anisotropy in the cosmic microwave background (CMB) have played 
a major role in establishing and sharpening the standard "ACDM" model of cosmology: a six-parameter model based 
on a fiat universe, dominated by a cosmological constant. A, and cold-dark-matter (CDM), with initial Gaussian, 
adiabatic fiuctuations seeded by inflation. This model continues to describe all existing CMB data, including the 
Wilkinson Microwav e Anisotropy Probe [WMAP) nine- year data presented in this paper and its companion paper 
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(pBennett et al. 2012), the small-scale temperature data ( Das et al. |2011 Keisler et al... 
Story et ai. ||2012 ), and the small-scale polarization data (Brown et al. 120091" Chiang et al. 12010 
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Reichardt et al. 2012a 
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Despite its notable success at describing all current cosmological data sets, the standard model raises many questions: 
what is the nature of dark matter and dark energy? What is the physics of inflation? Further, there are open questions 
about more immediate physical parameters: are there relativistic species present at the decoupling epoch, beyond the 
known photons and neutrinos? What is the mass of the neutrinos? Is the primordial helium abundance consistent 
with Big Bang nucleosynthesis? Are the initial fluctuations adiabatic? Tightening the limits on these parameters is 
as important as measuring the standard ones. Over the past decade WMAP has provided a wealth of cosmological 
information which can be used to address the above questions. In this paper, we present the final, nine-year constraints 
on cosmological parameters from WMAP. 
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The paper is organized as follows. In Section[2[ we briefly describe the nine- year WMAP likelihood code, the external 
data sets used to complement WMAP data, and we update our parameter estimation methodology. Section |3] presents 
nine-year constraints on the minimal six-parameter ACDM model. Section[4]presents constraints on parameters beyond 
the standard model, such as the tensor-to-scalar ratio, the running spectral index, the amplitude of isocurvature modes, 
the number of relativistic species, the mass of neutrinos, spatial curvature, the equation-of-state parameters of dark 
energy, and cosmological birefringence. In Section [Sj we discuss constraints on the amplitude of matter fluctuations, 
CTg, derived from other astrophysical data sets. Section [6] compares WMAP^s seven-year measurements of the Sunyaev- 
Zel'dovich effect with recent measurements by Planck. In Section [Tj we update our analysis of polarization patterns 
around temperature extrema, and we conclude in Section |8] 

2. METHODOLOGY UPDATE 

Before discussing cosmological parameter fits in the remainde r of the paper, we summarize changes in our parameter 
esti mat ion methodology and our choice of input data sets. In §2.1| we review changes to the WMAP likelihood code. 
In j ]2.2| we discuss our choice of external data sets used to cornplement WMAP data in various tests. Most of these 
data sets are new since the seven-year data release. We conclude with some updates on our implementation of Markov 
Chains. 

2.1. WMAP Likelihood Code 

For the most part, the structure of the likelihood code remains as it was in the seven-year WMAP data release. 
Howe ver, i nstead of using the Monte Carlo Apodiscd Spherical Transform EstimatoR (MASTER) estimate (Hivon] 
et al.|[2002] ) for the I > 32 TT spectrum, w e now use an optimally-e stimated power spectr um and errors based on the 



quadratic estimator from Tegmark (19971, as discussed in detail in Bennett et al. (2012 1. This I > 32 TT sp ectrum 



is based on the template-cleaned V- and W-band data, and the KQ85y9 sky mask (see iBennett et al. (|2012|) for an 
update on the analysis ma sks). The likeliho od fu nction for / > 32 continues to use the Gaussian plus log-normal 
approximation described in Bond et al. (1998|) and Verde et al. (20031. 

The / < 32 TT spectrum uses the Blackwefl-Rao estimator, as before. This is bas ed on Gibbs samples obtained from 
a nine- year one-region bias-corrected ILC map described in ([ Bennett et aL 2012) and sampled outside the KQ85y9 
sky mask. The map and mask were degraded to HEALPix \ and 2 /zK of random noise was added to each pixel in 
the map. 

The form of the polarization likelihood is unchanged. The I > 23 TE spectrum is based on a MASTER estimate 
and uses the template-cleaned Q-, V-, and W-band maps, evaluated outside the KQ85y9 temperature and polari zation 
masks. The / < 23 TE, EE, and BB likelihood retains the pixel-space form described in Appendix D of Page et al.| 
([2007). The inputs arc template-cleaned Ka-, Q-, and V-band maps and the HEALPix r3 polarization mask used 
previously. 

As previously, the likelihood code accounts for several important effects: mode coupling due to sky masking and 
non-uniform pixel weighting (due to non-uniform noise) ; beam window function uncertainty, which is correlated across 
the entire spectrum; and resid ual point source su b traction uncertaint y , which is also highly correlated. The treatment 
of these effects is described in [Verde et al.[ ([2003[); [Nolta et al.[ (|2009[); Punkley et al.[(2009[). 



2.2. External Data Sets 
2.2.1. Small-scale CMB measurements 

Since the time when the seven-year WMAP analyses were published, there have been new measurem ents of small- 
scale CMB fluctuations by th e Atacama Cosmology Telescope (ACT) ( Fowler et al. 2010 Das et al. ^2011) and the 



South Pole Telescope (SPT) (Keisler et al. 2011 Reichardt et al. 20l2ar 'I'hey have reported the angular power 
spectrum at 148 and 217 GHz for ACT, and at 95, 150, and 220 GHz for SPT, to 1' resolution, over ~1000 deg^ 
of sky. At least seven acoustic peaks are observed in the angular po wer spectrum, and the results are in remarkable 
agreement with the model predicted by the WMAP seven-year data ( Keisler et al.[[20lT ). 

Figure [l] shows data from ACT and SPT at 150 GHz, which constitutes the exte nded CMB data set used extensively 
in this paper (subsequently denoted 'eCMB'). We incorporate the SPT data from Keisler et al. (2011 1, using 47 band- 
powers in the range 600 < I < 3000. The likelihood is assumed to be Gaussian, and we use the published band-power 
window function s an d covariance matrix, the latter of which accounts for noise, beam, and calibration uncertainty. 
As discussed in j ]2.3[ we account for residual extragalactic foregrounds by marginalizing ove r three parameters: the 
Poisson and clustered point s ource amphtudes , and the SZ amplitude (Keisler et al. 2011). For ACT we use the 
148 GHz power spectrum from Das et al. (2011) in the multipole range 5(J0 < / < 10000, marginalizing over the same 
clustered point source and SZ amplitudes as m the SPT likelihood, but over a separate Poisson source amplitude. 

In ad dition to the temperature spectra, both AC T and SPT have estimated the deflection spectra due to gravitational 
lensing ( Das et aL[ 2011 van Engelen et aL[[2012 ). These measurements are consistent with predictions of the ACDM 
model fit to WmAp. Wlien we incorporate SP'l and ACT data in the nine-year analysis, we also include the lensing 
likelihoods provided by each groufFjto further constrain parameter fits. 

New observa.tions of the CMB polarization power spectra have also been released by the QUIET experiment ( QUIET 
Collaboration 2011 2012); their TE and EE polarization spectra are in excellent agreement with predictions based 



^® The map resolution levels refer to the HEALPix pixelization scheme ( jCorski et al.|2005| where r4, r5, r9, and rlO refer to A^aide values 
of 16, 32, 512, and 1024, respectively. ^ 
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these codes are available at http://lainbda.gsfc.nasa.gov 
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primarily on WMAP temperature fluctuation measurements. These data are the most recent in a series of polarization 
measurements at I > 50. However, high-? polarization observations do not (yet) substantially enhance the power of 
the full data to constrain parameters, so we do not include them in the nine-year analysis. 
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Fig. 1. — A compilation of the CMB data used in the nine-year WMAP analysis. The WMAP data are shown in black, the extended 
CMB data set — denoted 'eCMB' throughout — includes SPT data in blue {Keisler et al. 2011), and ACT data in orange, (Das et al. 20111. 
We also incorporate constraints from CMB lensing published by the SPT and AU'l' groups (not shown). The ACDM model tit to tne 
WMAP data alone (shown in grey) successfully predicts the higher-resolution data. 



2.2.2. Baryon Acoustic Oscillations 

The acoustic peak in the galaxy correlation function has now been detected over a range of redshifts from z — OA 
to z = 0.7. This linear feature in the galaxy data provides a standard ruler with which to measure the distance ratio, 
Dv /ts, the distance to objects at redshift z in units of the sound horizon at recombination, independent of the local 
Hubble constant. In particular, the observed angular and radial BAG scales at redshift z provide a geometric estimate 
of the effective distance, 

Dv{z) 



[{l + zfD\{z)czlH{z)Y/\ 



(1) 



where Da{z) is the angular diameter distance and H{z) is the Hubble parameter. The measured ratio Dy |rs^ where 
is the co-moving sound horizon scale at recombination, can be compared to theoretical predictions. 
Since the release of the seven-year WMAP data, the acoustic scale has been more precisely measured by the Sloan 
Digital Sky Survey (SDSS) and SDSS-HI Baryon Oscillation Spectroscopic Survey (BOSS) galaxy surveys, and by 
the WiggleZ and 6dFGS surveys. Previously, ov er half a million galax ies and LRGs from the SDSS-DR7 catalog had 

0.2 and 



been combined with galaxies from 2dFGRS by Percival et al 
0.35 



(20101 to measure the acoustic scale at 



et al 



(These data were used in the WMAP seven-year analysis.) Using the reconstruction method of [ Eisenstein 
;e of the acoustic scale in the SE 

± 0.17, and reducing the uncertainty from 3.5% to 1.9%. More recently the SUSS 



L] (2007), an improved estimate of the acoustic scale in the SDSS-DR7 data was made by [Padmanabhan et al 
pUT^ ^lrnhg Dy(0.35)/r 



L)K9 data from the BOSS survey has been used to estimate the BAO scale of the GMASS sample 
Dv{0 .^1) jr s = 13.67±0.22 for galaxies in the range 0.43 < z < 0.7 (at an eff e ctive redshift z 
20121. This result is used to constrain cosmological models in Sanchez et al. (2012). 



They report 



0.57) (Anderson et al. 
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TABLE 1 

BAO Data Used in the Nine- year. Analysis 



Redshift 



Data Set 



rs/Dv{z) 



Ref. 



0.1 6dFGRS 0.336 ±0.015 [Beutler et al| ( |2011^ 



0.35 
0.57 

0.44 
0.60 
0.73 



SDSS-DR7-rec 0.113 ± 0.002'" iP admanabhan et al. | || 2012[) 
SDSS-DR9-rec 0.073 ± 0.001=" [Anderson et aL| | |aJl2| | 



WiggleZ 
WiggleZ 
WiggleZ 



0.0916 ±0.0071 
0.0726 ± 0.0034 
0.0592 ± 0.0032 



P 



Blake et iiV^WV^ 
lake et al. (2012) 



(Blake 



et al., 1 ,2012 > 



^ For uniformity, the SDSS values given here have been inverted from the 
published values: Dv(0.35)/rs = 8.88±0.17, and Dv(0.57)/rs = 13.67±0.22. 



The acoustic scale has also been measured at higher redshift using the WiggleZ galaxy survey, 
report distances in three correlated redshift bins betwe en 0.44 and 0.73. A t lower redshift, 



Blake et al.| ( 2012 | 
0.1, a detection of 



the BAO scale has been made using the 6dFGS survey ( [Beu tler et al.ll2011|). These measurements are summarized in 
Table [Tj and plotted as a function of redshift in Figure 19 of Anderson ct al. (2012), together with the best-fit ACDM 



model prediction from the WMAP seven-year analysis ( [Komatsu et ai.„2011 
CMB-based prediction over the measured redshift range. 
For the nine- year analysis, we incorporate these data into a likelihood of the form 



I'he BAO data are consistent with the 



21nL = (x-d)^C-i(x-d) 



where 



and 



x-d=[r,/i^y(0.1)- 



0.336, Dv{<d.?,^)/r, - 
- 0.0916, rs/Dv{Q.m) - 



.88, L>y(0.57)/r^ - 13.67, 
0.0726, rslDv{Q.7i) - 0.0592] 



/ 4444.4 \ 

34.602 * 

20.661157 

24532.1 -25137.7 12099.1 

GO -25137.7 134598.4 -64783.9 

V 12099.1 -64783.9 128837.6 / 



(2) 



(3) 



(4) 



The mod el distances are derive d from the ACDM parameters using the same scheme we used in the WMAP seven-year 
analysis ( Komatsu et al.|[2011 ). 



2.2.3. Hubble Parameter 



It is instructive to combine WMAP measurements with measurements of the current expansion rate of the universe. 
Recent advances in the determinatio n of the Hubble constant have been made si nce the two teams using HST/WFPC2 



observations reported thei r results (Freedman et al. 2001 Sandage et al. 2006). Re-anchoring the HST Key Project 
distance ladder technique, Freedman et al. "( |2012[ ) report a significantly improved result of Hq = 74.3 ± 1.5 (statistical) 
±2.1 (systematic) km sec~^ Mpc~^. The overall 3.5% uncertainty must be taken with some caution however, since 
the uncertainties in all r ungs are not fully propagated. 

In a parallel approach, Riess et al. ( 2009 1 redesigned the distance ladder and its observations to control the systematic 
errors that dominated tne measurements. These steps include: the elimination of zero-point uncertainties by use of 
the same photometric system across the ladder; observations of Cepheids in the near-infrared to reduce extinction 
and sensitivity to differences in chemical abundance (the so-called "metallicity effect" ) ; the use of geometric distance 
measurements to provide a reliable absolute calibration; and the replacement of old Type la supernovae observations 
with recent ones that use the same photometric sys tems that define the Hubble flow. This approach has led to a 



measurement of Hq ~ 73.8 ± 2.4 km sec Mpc (Riess et al. 2011) with a total uncertainty of 3.3%. Since this 



uncertainty is smaller and has more fully propagated uncertainties, we adopt it in our analysis. 

2.2.4. Type la Supernovae 

The first direct evidence for acceleration in the expansion of the universe ca me from m easurements of luminosity 
distance as a function o f redshift using Type la supernovae as standard candles (Riess et al. 1998; Schmidt ct al._^1998j 



Perlmutter et al.] 1999 ) . Numerous follow-up observations have been made, extending these early measurements to 
higher redshitt. After the seven-year WMAP analysis was published, the Supernova Legacy Survey analyze d their 
three-year sample ('SNLS3') of high redshift supernovae (Guy et al. 2010 Conley et al. 2011 Sullivan et al. 2011 ). They 
measured 242 Typ e la supernovae in the redshift range L).08 < z < 1.06, three times more than their firs t-year sample 
([Astier et ar]|2006|). The SNLS team combined the 242 SNLS3 supernovae with 123 SNe at low redshift ([Hamuy et aT 
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1996| iRiess et al.||1999 



searcnl Holtzman et a 



Jha et al.||2006[ [Hicken et al.||2009 
2UU8| ), and 14 blNe at z > 1 trom ' 



Contreras et al.||2010|) 93 SNe fro m the SDSS supernovae 



S i' measurements by Riess et al. ( 2007 ) to form a sample 



of 472 SJNe. All ot these supernovae were re-analyzed using both the SALT2 and SiFTD light curve fitters, which give 
a n estimate of the SN peak rest-frame B-band apparent magnitude at the epoch of maximum light in that filter. 



Sullivan et al. (2011) carry out a cosmological analysis of this combined data, accounting for systematic uncertainties 
including common photometric zero-point errors and selection effects. They adopt a likelihood of the form 



21nL = (ms - m%YC^\^B ~ m 



(5) 



where is the peak- light apparent magnitude in B-band for each supernova, is the corresponding magnitude 
predicted by the model, and C is the covariance matrix of the data. Their analysis assigns three terms to the covariance 
matrix, C = Dgtat + Cstat + Csys, where Dgtat contains the independent (diagonal) statistical errors for each supernova, 
Cstat includes the statistical errors that are correlated by the light-curve fitting, and Cgys has eight terms to track 
systematic uncertainties, including calibration errors. Milky Way extinction, and redshift evolution. The theoretical 
magnitude for each supernova is modeled as 



51og[dL(z)/Mpc] - a(s -1)+I3C + M, 



(6) 



where Ai the empirical intercept of the — z relation. The parameters a and f3 quantify the stretch-luminosity and 
color-l uminosity relationship s, and the statistical error, Cgtat, is coupled to both parameters. Assuming a constant w 
model, Sullivan et al. (20111) measure a = 1.37 ± 0.09, and f3 — 3.2 ± 0.1. Including the term Cgys has a significant 
effect: it mcreases the error on the dark energy equation of state, from 0.05 to 0.08 in a flat universe. 

The 472 Type la su pernovae used in the SNLS3 analysis are consistent with the ACDM model predicted by WMAP 
( Sullivan et al.|[2011 ), thus we can justify including these data in the present analysis. However, the extensive study 
presented by the SJNLS team shows that a significant level of systematic error still exists in current supernova observa- 
tions. Hence we restrict our use of supernova data in this paper to the subset of models that examine the dark energy 
equation of state. When SNe data are included, we marginalize over the three parameters a, f3, and M. a and are 
sampled in the Markov Chain Monte Carlo (MCMC) chains, while M is marginalized analytically (Lewis & Bridle 
[20G2| ). 

2.3. Markov Chain Methodology 

As with previous WMAP analyses, we use MCMC methods to evaluate the likelihood of cosmological parameters. 
Aside from incorporating new likelihood codes for the external data sets described above, the main methodological 
update for the nine-year analysis centers on how we marginalize over SZ and point source amplitudes when analyzing 
multiple CMB data sets (i.e., " W MAP -|-eCMB"). W e have also incorp orated updates to the Code for Anisotropics in 
the Microwave Background (CAMB, 'L ewis et a"L]|2000 l, as described in S 2.3.1 

SZ amplitude - When combining data from multiple CMB experiments {WMAP, ACT, SPT) we sample and 
marginalize over a single SZ amplitude, Agz, that parameterizes the SZ contribution to all t hree data sets. To do so , 
we adopt a common SZ power spectrum template, and scale it to each experiment as follows. Battaglia et al. (2012b I 
compute a nominal SZ power spectrum at 150 GHz for the SPT experiment (their Figure 5, left panel, blue curve) . We 
adopt this curve as a spectral template and scale it by a factor of 1.05 and 3.6 to describe the relative SZ contribution 
at 148 GHz (for ACT) and 61 GHz (for WMAP), respectively. The nuisance parameter Asz then multiplies all three 
SZ spectra simultaneously. 

The above frequency scaling assumes a thermal SZ spectrum. For WMAP we assume an effective frequency of 61 
GHz, even though the WMAP power spectrum includes 94 GHz data. We ignore this error because WMAP data 
provide negligible constraints on the SZ amplitude when analyzed on their own. In the SPT and ACT frequency range, 
the thermal SZ spectrum is very similar to the kinetic SZ spectrum, so our procedure effectively accounts for that 
contribution as well. 

Clustered point sources - We adopt a common parameterization for the clustered point source contribution to 
both the ACT and SPT data, namely 1(1 + l)C//27r = ^cps l"'^ ( [Addison et"al]|2012[ ). Both the ACT and SPT teams 
use this form in their separate analyses at high I. (At low I, the SF'i' group adopts a constant spectrum, but this 
makes a negligible difference to our analysis.) By using a common amplitude for both experiments, we introduce one 
additional nuisance parameter. 

Poisson point sources - For unclustered residual point sources we adopt the standard power spectrum C; = const. 
Since the ACT and SPT groups use different algorithms for identifying and removing bright point sources, we allow the 
templates describing the residual power to have different amplitudes for the two experiments. This adds two additional 
nuisance parameters to our chains. 



2.3.1. CAMB 



Model po wer spectra are computed using the Code for Anisotropics in the Microwa ve Background (CAMB, Lewis 
et al. 2000), which is based on the earlier code CMBFAST (Seljak & Zaldarriaga 1996). We use the January 2012 
version ot CAMB throughout the nine- year analysis except when evaluating the [wq, Wa) model, where we adopted the 
October 2012 version. We adopt the default version of reef ast that is included with CAMB instead of other available 
options. As in the seven-year analysis, we fix the reionization width to be Az = 0.5. Since the WMAP likelihood 
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TABLE 2 

Maximum Likelihood ACDM Parameters'* 



Parameter 


Symbol 


WMAP data 


Combined data'' 


Fit ACDM parameters 






Physical baryon density 




0.02256 


0.02240 


Physical cold dark matter density 




0.1142 


0.1146 


Dark energy density {w = —1) 


Oa 


0.7185 


0.7181 


Curvature perturbations, kg = 0.002 Mpc~^ 




2.40 


2.43 


Scalar spectral index 


Us 


0.9710 


0.9646 


Reionization optical depth 


T 


0.0851 


0.0800 


Derived parameters 






Age of the universe (Gyr) 


to 


13.76 


13.75 


Hubble parameter, Ho = lOO/i km/s/Mpc 


Ho 


69.7 


69.7 


Density fluctuations @ 8h~^ Mpc 


o-g 


0.820 


0.817 


Baryon density /critical density 




0.0464 


0.0461 


Cold dark matter density /critical density 


He 


0.235 


0.236 


Redshift of matter-radiation equality 




3273 


3280 


Redshift of reionization 


•^rcion 


10.36 


9.97 



^ The maximum-likelihood ACDM parameters for use in simulations. Mean parameter values, 
with marginalized uncertainties, arc reported in Table [4] 
"Combined data" refers to WMAP-|-eCMB-|-BAO+ffo- 

code only incorporates low I BB data (with low sensitivity), we set the accurate_BB flag to FALSE and run the code 
with = 0. We set the high_accuracy_def ault flag to TRUE. When calhng the ACT and SPT likehhoods, we set 
k_eta_max_scalar — 15000 and l_max_scalar = 6000. The ACT likelihood extends to I = 10000, but foregrounds 
dominate beyond I = 6000. Except when exploring neutrino models, we adopt zero massive neutrinos and 3.04 massless 
neutrinos. The CMB temperature is set to 2.72548 K (Fixscn 2009). 

3. THE SIX-PARAMETER ACDM MODEL 
In this section we discuss the determination of the standard ACDM para mete rs, first using only the nine- year WMAP 



data, then, in turn, combined with the additional data sets discussed in ^2.2 
Carlo Markov Chain (MCMC) formalism used in previous analyses ([Sperge 



et al. 1 120071 iDunkley et al.|l2009^ 'Komatsu et al.||2009l iLarson et al.||2011 



Our analys i s employs the same Monte 
et al.||2003[ jVerde et al .1[2003l ISpergel 



Komatsu et al. 2011[ ). This formalism 



naturally produces parameter likelihoods that are marginalized over all other fit parameters in the model. Throughout 
this paper, we quote best-fit values as the mean of the marginalized likelihood, unless otherwise stated (e.g., mode or 
upper limits). Lower and upper error limits correspond to the 16% and 84% points in the marginalized cumulative 
distribution, unless otherwise stated. 

The six parameters of the basic ACDM model are: the physical baryon density, ilbh^; the physical cold dark matter 
density, flch^, the dark energy density, in units of the critical density, the amplitude of primordial scalar curvature 
perturbations, at fc = 0.002 Mpc~^; the power-law spectral index of primordial density (scalar) perturbations, 
n^; and the reionization optical depth, r. In this model, the Hubble constant, Hq = lOOh km/s/Mpc, is implicitly 
determined by the flatness constraint, fib + Clc + = 1. A handful of parameters in this model take assumed values 
that we further test in f|4j other parameters may be derived from the fit, a s in Table [2} Throug hout this paper we 



assume the initial fluctuations are adiabatic and Gaussian distributed (see Bennett et al. (20121 for limits on non- 



Gaussian fluctuations from the nine- year WMAP data) except in [ 4.2 where we allow the initial fluctuations to include 
an isocurvature component 



To assess WMAP data consistency, we begin with a comparison of the nine-year and seven-year results (Komatsu 
et al. 2011); we then study the ACDM constraints imposed by the nine- year WMAP data, in conjunction with the 



most recent external data sets available. 



3.1. Comparison With Seven-year Fits 

Table |3] gives the best-flt ACDM parameters (mean and standard deviation, marginalized over all other parameters) 
for selected nine-year and seven-year data combinations. In the case where only WMAP data are used, we evaluate 
parameters using both the C~^-weighted spectrum and the MASTER-based one. For the case where we include BAO 
and Hn p riors, we use only the C~"'^-weighted spectrum fo r the nine-year WMAP data, and we update the priors, as 
per { 2.2 The seven-year results are taken from Table 1 of Komatsu et al. (2011 1. 



3.1.1. WMAP Data Alone 

We first compare seven-year and nine-year results based on the MASTER spectra. Table [3] shows that the nine- 
year ACDM parameters are all within 0.5(t of each other, with flch^ having the largest difference. We note that the 
combination flmh^ + is approximately constant between the two models, reflecting the fact that this combination is 
well constrained by primary CMB fluctuations, whereas — is less so due to the geometric degeneracy. Turning 
to the -weighted spectrum, we note that the nine-year ACDM parameters based on this spectrum are all within 
~ 0.3a of the seven-year values. Thus we conclude that the nine-year model flts are consistent with the seven-year flt. 
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Next, we examine the consistency of the two ACDM model fits, derived from the two nine-year spectrum estimates. 
As seen in Table [3j the six parameters agree reasonably well, but we note that the estimates for Us differ by 0.75a, 
which we discuss oelow. To help visualize the fits, we plot both spectra (C~^-weighted and MASTER), and both 
models in Figure [2j As noted in Bennett et al. (2012), the difference between the two spectrum estimates is most 
noticeable in the range / '-^ 30 — 60 where the C~ ^-weighted spectrum is lower than the MASTER spectrum, by up 
to 5% in one bin. However, the ACDM model fits only differ noticeably for / < 10 where the fit is relatively weakly 
constrained due to cosmic variance. 

To understand why these two model spectra are so similar, we examine parameter degeneracies between the six 
ACDM parameters when fit to the nine-year WMAP data. In Figure [3] we show the two largest degeneracies that 
affect the spectral index Ug, namely 10^ and Q,bh? . The contours show the 68% and 95% CL regions for the fits 
to the C~^-weighted spectrum while the plus signs show the maximum likelihood points from the MASTER fit. Note 
that the C~^-weighted fits favor lower 10^ and higher fib/i^, both of which push the C~^-weighted fit towards 
higher Ug. Given the consistency of the fit model spectra, we conclude that the underlying data are quite robust and 
in subsequent subsections, we look to external data to help break any degeneracies that remain in the nine-year data. 

We conclude this subsection with a summary of some additional tests with simulations that we carried out to assess 
the robustness of the C~^-weighted spectrum estimate in general, and the fits in particular. First, we compared 
C~^-weighted and MASTER spec tra obtained from the 500 "parameter recovery" simulations developed for our seven- 



year analysis (Larson et al. 20111. In this test, the C -weighted spectrum was shown to be unbiased over the full 



multipole range that WMAy is sensitive to 

We next evaluated a number of difference statistics, but the one that was deemed most pertinent to understanding 
the Us fit was the average power difference between / = 32 — 64 (this is admittedly a posterior choice of I range). 
When the parameter recovery simulations were analyzed with the conservative KQ75y9 mask, more than one-third of 
the simulated spectra had a larger power difference (C~^— MASTER, in the Z = 32 — 64 bin) than did the flight data. 
However, when the same analysis was performed with the KQ85y9 mask, only 2 out of 500 simulation realizations had 
a larger difference than did the flight data. This led us to investigate residual foreground contamination. 

We amended the CMB-only parameter recovery simulations with model foreground signals that we deemed to be 
representative of both the raw foreground signal outside the KQ85y9 mask, and an estimate of the residual contamina- 
tion after template cleaning. These signals were based on the modeling studies described in Bennett et al.| ([2012). We 
repeated the comparison of the two spectrum estimates with the foreground-contaminated simulations. Both estimates 
showed slightly elevated power in the Z = 32 — 64 bin (a few percent) , with the MASTER estimate being slightly more 
sensitive to residual foreground power. However, the distribution of spectrum differences was not significantly different 
than with the CMB-only simulations. In the end, we attribute the spectrum differences to statistical fiuctuations and 
we adopt the C^^-weighted spectrum for our final analysis because it has lower uncertainties (Bennett et al. 2012). 

To conclude the seven-year/nine- year comparison, wc note that the remaining 5 ACDM parameters changed by less 
than 0.3(7 indicating very good consistency. The overall effect of the nine-year WMAP data is to improve the average 
parameter uncertainty by about 10%, with Viji? and J^a each improvin g by nearly 20%. Th e latter improvement 
is a result of higher precision in the third acoustic peak measurement ( Bennett et al.| 2012| which gives a better 
determination of Vlch?. This, in t urn, improves JIa, which is constrained by flatness (or in non-flat models, by the 
geometric degeneracy discussed in ^4.5). The overall volume reduction in the allowed 6-dimensional ACDM parameter 



space in the switch from seven-year to nine-year data is a factor of 2, the majority of which derives from switching to 
the C~ ^-weighted spectrum estimate. 



3.1.2. WMAP Data With BAO and Hi 







To complete our comp arison with seven-year results, we examine ACDM fits that include the BAO and Hq priors. 
In Komatsu et al. (2011) we argued that these two priors (then based on earlier data) provided the most robust and 
complementary parameter constraints, when used to supplement WMAP data. In Table [3] we give results for the 
updated version of this data c ombin atio n, whic h includes the nine-year C~^-weighted spectrum for WMAP and the 
BAO and Hp priors noted in j|2.2.2| and |2.2.3[ respectively. For comparison, we reproduce seven-year numbers from 



Table 1 of [Komatsu et al.| ( [20Tl| ' 



As a measure of data consistency, we note that 4 of the 6 ACDM parameters changed by less than 0.25(7 (in units 
of the seven-year a) except for ftch^ and fl\ which changed by ±0.8(7, respectively. As noted above, these latter 
two parameters were more stable when fit to WMAP data alone (both in absolute value and in units of cr), so we 
conclude that this small change is primarily driven by the updated BAO and Hq priors. In particular, CMB data 
provide relatively weak constraints along the geometric degeneracy line (which corresponds to a line of nearly constant 
fljn + when spatial curvature is allowed), so external data are able to force limited anti-correlated changes in 
{n,n,n\) with relatively little penalty from the WMAP likelihood. In subsequent sections we explore the nine-year 
ACDM fits more fully by adding external data sets to the WMAP data one at a time. 

The combined effect of the nine-year WMAP data and updated the BAO and Hq priors is to improve the average 
parameter uncertainty by nearly 25%, with flch^ and f^A each improving by 37%, due, in part, to improved constraints 
along the geometric degeneracy line. The overall volume reduction in the allowed 6-dimensional ACDM parameter 
space is a factor of 5, nearly half of which (a factor of 2) comes from the nine-year WMAP data alone. 



3.2. ACDM Constraints From CMB Data 
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Fig. 2. — Two estimates of the WMAP nine-year power spectrum along with the best-fit model spectra obtained from each; black - the 
C~^-weighted spectrum and best fit model; red - the same for the MASTER spectrum and model. The two spectrum estimates differ by 
up to 5% in the vicinity of i ~ 50 which mostly affects the determination of the spectral index, rts, as shown in Table [s] We adopt the 
^-weighted spectrum throughout the remainder of this paper. 
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109A| Qbh2 

Fig. 3. — 68% and 95% CL regions for the ACDM parameters n^, 10^A|j, and Q,t,h? . There is a modest degeneracy between these 
three parameters in the six-parameter ACDM model, when fit to the nine-year WMAP data. The contours are derived from fits to the 
C^'^-weighted power spectrum, while the plus signs indicate the maximum likelihood point for the fit to the MASTER power spectrum. 
As shown in Figure |2] the two model produce nearly identical spectra. 
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TABLE 3 

WMAP Seven-year to Nine-year Comparison of the Six-Parameter ACDM Model^ 



Parameter 


Nine-year 


WMAP-on\y° 
Nine-year (MASTER)<= 


Seven-year 


WMAP+BAO+Ho^ 
Nine-year Seven-year 


Fit parameters 














0.02264 ± 0.00050 


0.02243 ± 0.00055 


02249+°"°°^^ 


0.02266 ± 0.00043 


0.02255 ± 0.00054 




0.1138 ±0.0045 


0.1147 ±0.0051 


0.1120 ±0.0056 


0.1157 ± 0.0023 


0.1126 ±0.0036 


Ha 

ris 

T 


0.721 ±0.025 

2.41 ±0.10 
0.972 ±0.013 
0.089 ±0.014 


0.716 ±0.028 

2.47 ±0.11 
0.962 ±0.014 
0.087 ± 0.014 


n -97+0.030 
'J- '^'-0.029 
2.43 ±0.11 

0.967 ± 0.014 
0.088 ±0.015 


0.712 ±0.010 

9 ,97+0.078 
^•^^'_0.079 

0.971 ± 0.010 
0.088 ±0.013 


0.725 ±0.016 
2.430 ±0.091 
0.968 ±0.012 
0.088 ±0.014 



Derived parameters 

to (Gyr) 

Ho (km/s/Mpc) 



13.74 ±0.11 
70.0 ±2.2 
0.821 ± 0.023 
0.0463 ± 0.0024 
0.233 ±0.023 
10.6 ± 1.1 



13.75 ± 0.12 
69.7 ±2.4 
0.818 ±0.026 
0.0462 ± 0.0026 
0.237 ± 0.026 
10.5 ± 1.1 



13.77 ± 0.13 
70.4 ±2.5 

r, oi 1+0.030 

U.«ii_o.031 

0.0455 ± 0.0028 
0.228 ±0.027 
10.6 ± 1.2 



13.750 ±0.085 
69.33 ± 0.88 
0.830 ±0.018 
0.0472 ± 0.0010 
0.2408_(, (jog2 
10.5 ± 1.1 



13.76 ±0.11 
70.2 ± 1.4 
0.816 ±0.024 
0.0458 ± 0.0016 
0.229 ±0.015 
10.6 ± 1.2 



^ Comparison of 6-parameter ACDM fits with seven-year and nine-year WMAP data, with and without BAO and Ho priors. 
^ The first three data columns giv e res ults from fitting to WMAP data only. The last two columns give results when BAO and Ho 
priors are added. As di scussed in ^^2?2\ these priors have been updated for the nine-year analysis. The seven-year results are taken 
directly from Table 1 of |Komatsu et a l. (2011). 

Unless otherwise noted, the nine-year WMAP likelihood uses the C~^-weighted power spectrum whereas the seven-year likelihood 
used the MASTER-based spectrum. The column labeled 'Nine-year (MASTER)' is a special run for comparing to the seven-year 
results. 

From the standpoint of astrophysics, primary CMB fluctuations, combined with CMB lensing, arguably provide the 
cleanest probe of cosmology because the fluctuations dominate Galactic foreground emission over most of the sky, 
and they can (so far) be understood in terms of linear perturbation theory and Gaussian statistics. Thus we next 
consider parameter constraints that can be obtained when adding addit ional CMB data to the nine- year WMAP data. 
Specifically, we examine the effects of adding SPT and ACT data (see ^ 2.2.11: the best- fit parameters are given in the 
"-l-eCMB" column of Table H 

With the addition of the high-^ CMB data, the constraints on the energy density parameters flbh^, ^Ji^ , and 17a 



all improve by 25% over the precision from WMAP data alone. The improvement in the baryon density measurement 
is due to more precise measurements of the Silk damping tail in the power spectrum at / > 1000; the improvements in 
^Ji^ and r^A are du e in part to improvements in the high-Z T T data, but also to the detection of CMB lensing in the 
SPT and ACT data ( |Das et al.|20lT [van Engelen et aL||2012 [), which helps to constrain by fixing the growth rate 
of structure between z = 1100 and z = 1 — 2 (the peak m the lensing kernel). Taken together, CMB data available at 
the end of the WMAP mission produce a 1.6% measurement of Q,i,h? and a 3.0% measurement of Vt^h? . 

The increased fc-space lever arm provided by the high-/ CMB data improves the uncertainty on the scalar spectral 
index by 25%, giving Ug — 0.9646 ± 0.0098, which implies a non-zero tilt in the primordial spectrum (i.e., < 1) at 



3.6a. We examine the implications of this measurement for inflation models in S4.1 

If we assume a flat universe, which breaks the CMB's geometric degeneracy, then CMB data alone now provide a 
2.3% measurem ent o f the Hubble parameter, Hq — 70.5 ± 1.6 km/s/Mpc, independent of the cosmic distance ladder. 
As discussed in j ]3.4[ this is consisten t with the recent d etermination of the Hubble parameter using t he co smic distance 
ladder: Hq = 73. 8± 2.4 k m/s/ Mpc ( Riess et al. [2011 1; we explore the effect of adding this prior in S3. 4 We relax the 
assumption of flatness in §4.5[ 

3.3. Adding BAO Data 

Acoustic structure in the large scale distribution of galaxies is manifest on a co-moving scale of 152 Mpc, where the 
evolution of matter fluctuations is largely within the linear regime. A numb er of authors hav e studied the degree to 
which the acoustic structure could b e perturbed by nonline ar evolution ( e.g., Seo fc Eisensteint2 005 ; Seo & Eiscnstcin 



2007; Jeong & Komatsu 2006, 2009^: 'Cro cce fc Scoccimair o 2008; Matsubara 20Wl'aru ya Hira matsu 2008; PacF 
manabhan &: White 2009J , and the effects are well below the current measurement uncertainties. Because it is base J 
on the same well-understood physics that governs the CMB anisotropy, we consider measurements of the BAO scale 
to be the next-most robust cosmological probe after CMB fluctuations. The ACDM parameters fit to CMB and BAO 
data are given in the "+eCMB+BAO" column of Table J4] 

Measurements of the tangential and radial BAO scale at redshift z measure the effective distance Dv{z), given 
in equation in units of the sound horizon rs{zd). This quantity is primarily sensitive to the total matter and 
dark energy oensities, and to the current Hubble parameter. Since the BAO scale is relatively insensitive to the 
baryon density, flbh^, this parameter does not improve significantly with the addition of the BAO prior. However, the 
low-redshift distance information imposes complementary constraints on the matter density and Hubble parameter, 
improving the precision on D.ch'^ from 3.0% to 1.6%, and on Hq from 2.3% to 1.2%. In the context of standard ACDM 
these improvements lead to a measurement of the age of the universe with 0.4% precision: to = 13.800 ± 0.061 Gyr. 
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TABLE 4 

Six-Parameter, ACDM Fit; WMAP plus External Data^ 



Parameter 



WMAP 



heCMB 



+eCMB+BAO 



(-cCMB+Ho 



+eCMB+BAO+Ho 



Fit parameters 



Derived parameters 

to (Gyr) 

Ho (km/s/Mpc) 

0-8 



^rcion 



0.02264 ± 0.00050 0.02229 ± 0.00037 0.02211 ± 0.00034 0.02244 ± 0.00035 



0.1138 ±0.0045 
0.721 ± 0.025 

2.41 ± 0.10 
0.972 ±0.013 
0.089 ±0.014 

13.74 ±0.11 
70.0 ±2.2 
0.821 ± 0.023 
0.0463 ± 0.0024 
0.233 ±0.023 
3265liO« 
10.6 ± 1.1 



0.1126 ±0.0035 
0.728 ±0.019 
2.430 ±0.084 

0.9646 ± 0.0098 
0.084 ±0.013 

13.742 ± 0.077 
70.5 ± 1.6 

0.810 ±0.017 
0.0449 ± 0.0018 

0.227 ±0.017 
3230 ± 81 
10.3 ± 1.1 



0.1162 ± 0.0020 
0.707 ±0.010 

2.484_Q 

0.9579lH°ol^ 
0.079+oo;i 



13.800 ± 0.061 
68.76 ±0.84 
r, 099+0.013 

0.04678 ± 0.00098 
0.2460 ± 0.0094 
3312 ±48 
10.0 ± 1.0 



0.1106 ±0.0030 
0.740 ±0.015 

9 sqfi+0.079 
z.oao_o Q7g 

q690+° °°^l 
u.aDyu_Q Q090 

0.087 ±0.013 



13.702 ± 0.069 
71.6 ± 1.4 

0.803 ±0.016 
0.0438 ±0.0015 

0.216 ±0.014 
3184 ± 70 
10.5 ± 1.1 



0.02223 ± 0.00033 
0.1153 ±0.0019 

n 7lQc:+0.0095 

u. / loj-o.ooge 
2.464 ± 0.072 
0.9608 ± 0.0080 
0.081 ± 0.012 



13.772 ± 0.059 
69.32 ± 0.80 

0.820j:Hl! 
0.04628 ± 0.00093 

U.^4U^_0 0087 

3293 ± 47 
10.1 ± 1.0 



ACDM model fit to WMAP nine-year data combined with a progression of external data sets. 



The addition of the BAO prior helps to break some residual degeneracy between the primordial spectral index, n^, 
on the one hand, and ftch^ and Hq on the other. Figure [i] shows the 2-dimensional parameter likelihoods for {ns,ilch'^) 
and (ns,i/o) for the three data combinations considered to this point. With only CMB data (black and blue contours) 
there remains a weak degeneracy between Ug and the other two. When the BAO prior is added (red), it pushes rich^ 
towards the upper end of the range allowed by the CMB, and vice versa for Hq. Both of these results push Ug towards 
the lower end of its CMB-allowable range; consequently, with the BAO prior included, the marginalized measurement 
of the primordial spectral index is Hs = 0.9579^g Qgg2 which constitutes a 5a measurement of tilt (n^ < 1) in the 



primordial spectrum. We discuss the implications of this measurement for inflation models in ^4.1 
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Fig. 4. — Measurements of the scalar spectral index with CMB and BAO data. Left to right - contours of (Dv(0.57)/rs,ns), {Ho,ns), 
{flch^ ,na). Black contours show constraints using WMAP nine-year data alone; blue contours include SPT and ACT data ( WMAP+eCMB); 
red contours add the BAO prior( WMAP+cCMB+BAO). The BAO prior provides an independent measurement of the low-redshift distance, 
Dy {z)/rs, which maps to constraints on flch^ and Ho ■ When combined with CMB data, the joint constraints require a tilt in the primordial 
spectral index {ns< 1) at the 5cr level. 



3.4. Adding Hq Data 

Measurements of the Hubble parameter using the cosmic distance ladder have a long history, and are subject to 
a variety of different systematic errors that have been steadily reduced over time. However, an accurate, direct 
measurement of the current expansion rate is vital for testing the validity of the ACDM model because the value 
derived from the CMB and BAO data is model-dependent. Measurements of Hq provide an excellent complement to 
CMB and BAO measurements. The Hq prior considered here has a precision that approaches the ACDM-based value 
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given above. Consequently, we next consider the addition of the Hg prior discussed in |2.2.3[ without the inclusion of 
the BAO prior. The ACDM parameters fit to CMB and Hq data are given in "+eCMB+Bo column of Table |4] 

Two cosmological quantities that significantly shape the observed CMB spectrum are the epoch of matter radiation 
equality, z^q, which depends on Q,ch?, and the angular diameter dista nce to the last scattering surface, which 
depends primarily on Hq. As illustrated in Figure [5] (see also ^ 4.3.1), the CMB data still admit a weak degeneracy 
between ^ch^ and that the BAO and Hq priors help to break. The black contours in Figure [s] show the constraints 
from CMB data (WMAP+eCMB), the red from CMB and BAO data, and the blue from CMB with the Hq prior. 
While these measurements are all consistent, it is interesting to note that the BAO and Hq priors are pushing towards 
opposite ends of the range allowed by the CMB data for this pair of parameters. Given this minor tension, it is worth 
examining independent sets of constraints that do not share common CMB data. A simple test is to compare the 
marginalized constraints on the Hubble parameter from the CMB+BAO data {Hq — 68.76 ± 0.84 km/s/Mpc), to the 
direct, and independent, measurement from the distance ladder [Hq = 73.8 ± 2.4 km/s/Mpc). In our Markov Chain 
that samples the ACDM model with the WMAP+eCMB+BAO data, we found that only 0.1% of the values in 
the chain fell within the la range of the Hubble prior, but that 45% fell within the 2a range of 73.8 ± 4.8 km/s/Mpc. 
Based on this, we conclude that these measurements do not disagree, and that they may be combined to form more 
stringent constraints on the ACDM parameters. 

We conclude this subsection by noting that measurements of the remaining ACDM parameters are modestly improved 
by the addition of the Hq prior to the CMB data, with Vlch^ improving the most due to the effect discussed above and 
illustrated in Figure [5] 
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Fig. 5. — Measurements of Q,ch? and Hq from CMB data only (blue contours, WMAP+eCMB), from CMB and BAO data (green 
contours, WMAP+eCMB+BAO), and from CMB and Hq data (red contours, WMAP+eCMB+_ffo)- The two non-CMB priors push the 
constraints towards opposite ends of the range allowed by the CMB alone, but they are not inconsistent. 



3.5. ACDM Fits to the Combined Data 

Given the consistency of the data sets considered above, we conclude with a summary of the ACDM fits derived 
from the union of these data. The marginalized results are given in the "+eCMB+BAO+i/o" column of Table |4j 
The matter and energy densities, flbh^, ^ch^, and ft a are all now determined to ^1.5% precision with the current 
data. The amplitude of the primordial spectrum is measured to within 3%, and there is now evidence for tilt in the 
primordial spectrum at the 5a level. 

At the end of the WMAP mission, the nine-year data produced a factor of 68,000 decease in the allowable volume 
of the six-dimensional ACDM parameter space, relative to the pre- WMAP measurements. When these data are 
combined with the eCMB+BAO+i?o priors, we obtain an additional factor of 27 over the WMAP-only constraints. 
As an illustration of the predictive power of the current data. Figure [6] shows the la range of high-1 power spectra 
allowed by the six-parameter fits to the nine-year WMAP data. As shown in Figure [l] this model has already predicted 
the current small-scale measurements. If future measurements of the spectrum, for example by Planck, lie outside this 
range, then either there is a problem with the six-parameter model, or a problem with the data. 

Remarkably, despite this dr amatic increase in pre cision, the six-parameter ACDM model still produces an acceptable 



fit to all of the available data. Bennett et al. 
WMAP data. In the next section we place 



(2012) present a detailed breakdown of the goodness of fit to the nine-year 



imits on parameters beyond the six required to describe our universe. 



4. BEYOND SIX-PARAMETER ACDM 
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Fig. 6. — The nine-year WMAP data (in black) are shown with the Icr locus of six-parameter ACDM models allowed by the nine-year 
WMAP data. The error band is derived from the Markov Chain of six-parameter model fits to the WMAP data alone. The blue curve 
indicates the mean of the ACDM model fit to the WMAP-|-eCMB data combination. At high-/ this curve sits about la below the model 
fit to WMAP data alone. The marginalized parameter constraints that define these models are given in the WMAP and WMAP+eCMB 
columns of Table |4] 



In this section we discuss constraints that can be placed on cosmological param eters beyond the standard model 
using the nine-year WMAP data combined with the external data sets discussed in j ]2.2| In the following subsections, 
we consider limits that can be placed on additional parameters one or two at a time, beginning with constraints on 
initial conditions and proceeding through to the late-time effects of dark energy. 



4.1. Primordial Spectrum and Gravitational Waves 

As noted in §3.5[ the nine-year WMAP data, when combined with eCMB, BAO and Hq priors, exclude a scale- 
invariant primordial power spectrum at 5(t significance. For a power-law spectrum of primordial curvature perturba- 
tions, 

,{k) = A'^{ko)(-) , (7) 



with ko = 0.002 Mpc~^, we find = 0.9608 ± 0.0080. This result assumes that tensor modes (gravitational waves) 
contribute insignificantly to the CMB anisotropy. 

At this time, the most sensitive limits on tensor modes are still obtained from the shape of the temperature power 
spectrum, in conjunction with additional data. For example, Story et al. (2012) report r < 0.18 (95% CL), where 



A?(fco) 



_ Phjko) 
Aliko) Pniko)- 



(8) 



Due to confusion f rom density fluctuations, the lowest tensor amplitude that can be reliably detected from temperature 
data is r < 0.13 (Knox 1995). Several recent ex periments are b eginning to establish comparable limits from non- 
detection of B-mode polarization anisotropy, e.g., Chiang et al. ( [2010| ) report r < 0.7 (95% CL) from BICEP and 
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the QUIET Collaboratioii] ( 2012 1 reports r < 2.8 (95% CL). A host of forthco ming experiments are targeting B-mode 
measurements that have the potential to detect or limit tensor modes at significantly lower levels than can be achieved 
with temperature data alone. 

In Table [5j we report limits on r from the nine-year WMAP data, analyzed alone and jointly with external data; 
the tightest constraint is 

r < 0.13 (95% CL) WMAP+eCMB+BAO+Ho- 

This is effectively at the limit one can reach without B-mode polarization measurements. The joint constraints on 
and r are shown in Figure [7j along with selected model predictions derived from single-field inflation models. Taken 
together, the current data strongly disfavor a pure Harrison-Zel'dovich (HZ) spectrum, even if tensor modes are allowed 
in the model fits. 
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Fig. 7. — Two-dimensional marginalized constraints (68% and 95% CL) on the primordial tilt, Us, and the tensor-to-scalar ratio, r, derived 
with the nine-year WMAP in conjunction with: eCMB (green) and e CMB-|-BAO-| -f/n (red). The symbols and lines show predictions from 
single-field in flation models whose p otential is given by V(4>) oc (f)" ( |Linde||l983| l, with q = 4 (solid), a = 2 (long-dashed), and a = 1 
(short-dashed; [McAllister et al. |2010[|. Also shown are those from the first inflation model, which is based on an term in the gravitational 
Lagrangian (dotted; |Starobins^^^80| . Starobinsky's model gives ns = 1 — 2/N and r = 12/N'^ where N is the number of e-folds between 
the end of inflation and the epoch at which the scale k = 0.002 Mpc~^ left the horizon during inflatio n. These pre dictions are the same 
as those of inflation models with a £,4'^R term in the gravitational Lagrangian with a potential (Komatsu & Futamasc 1999). See 
Appendix \K\ for details. 



4.1.1. Running Spectral Index 

Some inflation models predict a scale dependence or "running" in the (nearly) power-law spectrum of scalar pertur- 
bations. This is conveniently p a.rameterized by the logarit hmic derivative of the spectral index, dns/dlnk, which gives 
rise to a spectrum of the form ( Kosowsky fc Turner||1995 ) 



/ U \ ^^s ('i^o)~l+ ^ ln(fc/fco)t^ns /rf In 

Alik) ^ Aliko) [^-j . (9) 

We do not detect a statistically significant deviation from a pure power-law spectrum with the nine-year WMAP data. 
The allowed range of dns/dlnk is both closer to zero and has a smaller confidence range with the nine- year data, 
dus/dlnk = —0.019 ± 0.025. However, with the inclusion of the high-/ CMB data, the full CMB data prefer a slightly 
more negative value, with a smaller uncertainty, dus/dlnk — — 0.022j^Q Q![^. While not significant, this result might 
indicate a trend as the Z-range of the data expand. The inclusion of BAO and Ho data does not affect these results. 

If we allow both tensors and running as additional primordial degrees of freedom, the data prefer a slight negative 
running, but still at less than 3cr significance, and only with the inclusion of the high-/ CMB data. Complete results 
are given in Table [5] 

4.2. Isocurvature Modes 

In addition to adiabatic fluctuations, where all species fluctuate in phase and therefore produce curvature fluctuations, 
it is possible to have isocurvature perturbations: an over-density in one species compensates for an under-density in 
another, producing no net curvature. These entropy, or isocurvature perturbations have a measurable effect on the 
CMB by shifting the acoustic peaks in the power spectrum. For cold dark matter and photons, we define the entropy 
perturbation field 
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TABLE 5 

Primordial spectrum: tensors & running scalar index^ 



Parameter WMAP +cCMB +eCMB+BAO +eCMB+BAO+i?() 

Tensor mode amplitude*^ 

< 0.13 (95% CL) 
0.9636 ± 0.0084 

-0.023 ±0.011 
1.020 ± 0.029 

< 0.47 (95% CL) 
-0.040 ±0.016 

1.075 ± 0.046 

^ A complete list of parameter values for these models may be found at 
http : //lambda . gsf c . nasa . gov/ 

The tensor mode amplitude and scalar running index parameter are each fit singly, and then 
jointly. In models with running, the nominal scalar index is quoted at ko = 0.002 Mpc~^. 



r < 0.38 (95% CL) < 0.17 (95% CL) < 0.12 (95% CL) 

Us 0.992 ±0.019 0.970 ±0.011 0.9606 ± 0.0084 

Running scalar index'' 

dris/dlnk -0.019 ± 0.025 -0.022 -0.024 ±0.011 

Us 1.009 ±0.049 1.018 ±0.029 1.020 ±0.029 

Tensors and running, jointly'' 

r < 0.50 (95% CL) < 0.53 (95% CL) < 0.43 (95% CL) 

dnjdlnk -0.032 ± 0.028 -0.039 ± 0.016 -0.039 ± 0.015 

ns 1.058 ±0.063 1.076 ±0.048 l.Oestnnli 



TABLE 6 

ISOCURVATURE MODES'' 



Parameter 



WMAP 



+cCMB 



t-eCMB+BAO +eCMB+BAO+Ho 



Anti-correlated modes'' 

a_i < 0.012 (95% CL) 

Uch"^ 0.1088 ±0.0050 

Us 0.994 ±0.017 

n 007+0.025 

Uncorrelated modes'^ 

ao < 0.15 (95% CL) 

Qch"^ 0.1093 ±0.0056 

Us 0.994 ±0.021 

as 0.805 ± 0.027 



< 0.0076 (95% CL) < 0.0035 (95% CL) 



0.1097 ± 0.0037 
0.977 ± 0.011 
0.802 ±0.018 

< 0.061 (95% CL) 
0.1115 ±0.0036 
0.970 ± 0.011 
0.805 ±0.018 



0.1160 ± 0.0020 
9631+°"°*'' 

r, 000 + 0.014 

U.»^3_0.oi3 



< 0.043 (95% CL) 
0.1161 ± 0.0020 
0.9608"* 



j+0.0086 
-0.0085 

0.821 ± 0.014 



< 0.0039 (95% CL) 
0.1151 ± 0.0019 

9662+° °°*'5 

n QOI +0.014 
U.»^l_0.013 

< 0.047 (95% CL) 
0.1152 ±0.0019 

9639+° °°*'' 
0.819 ±0.014 



A complete lis t of parameter values for these models may be found 



at 



Ihttp : //lambda . gsf c . nasa . gov/ 

The anti-correlated isocurvature amplitude comprises one additional parameter in the ACDM fit. 
The remaining parameters in this table section are given for trending. 

The uncorrelated isocurvature amplitude comprises one additional parameter in the ACDM fit. The 
remaining parameters in this table section are given for trending. 



( Bean et al.||2006 Komatsu et al.|[2009 |. The relative amplitude of its power spectrum is parameterized by a, 

a ^ Psjko) 
l-a~ Pniko)' 



(11) 



with fco = 0.002 Mpc-^ 

We consider two types of isocurvature modes: those which are completely uncorrelated with the curvature modes 
(with amplitude ao), motivated by the axion model, and those which are anti-correlated with the curvature modes (with 
amplitude a_i), motivated by the cur vaton model. For the latter, we adopt the convention in which anti-correlation 
increases the power at low multipoles ( Komatsu et al.|[2009 ). 

The constraints on both types of isocurvature modes are given in Table [6l We do not detect a significant contribution 
from either type of perturbation in the nine- year data, whether or not addition al data are include d in the fit. With 
WMAP data alone, the limits are slightly improved over the seven-year results (|Larson et al."20 11'), but the addition 
of the new eCMB data improves limits by a further factor of ~2. Adding the BAU data [\ 2.2.2 ) and Hq data ( p. 2.3 1 
further improves the limits, to 



a_i < 0.0039 (95% CL) 
ao < 0.047 (95% CL) 



WMAP+eCMB+BAO+Ho, 



due to the fact that these data help to break a modest degeneracy in the CMB anisotropy between the isocurvature 
modes and the ACDM parameters given in Table [6j 
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4.3. Number of Relativistic Species 
4.3.1. The number of relativistic species and the CMB power spectrum 
Let us write the energy density of relativistic particles near the epoch of photon decoupling, z ; 



1090, as 



P-1 + Pu + Per, 



(12) 



where, in natural units, = j^T^ is the photon energy density, pi, — -^j^N^T^ is the neutrino energy density, and 
Per denotes the energy density of "extra radiation species." (The factor of 7/8 in t he neutrino density arises from the 

In the standard model of particle physics 



Fermi-Dirac distribution. 



3.04 (Dicus et al. 



1972). 



1982 Mangano et al. 



2002 1 , while in the standard thermal history of the universe, = (4/11)1/3 (e.g., |Weinberg| 

Since we don't know the nature of an extra radiation species, we cannot specify its energy density or temperature 
uniquely. For example, por could be comprised of bosons or fermions. Nevertheless, it is customary to parameterize 
the number of extra radiation species as if they were neutrinos, and write 



Pi^ 



7^ 
120' 



p4 



(13) 



where iVcft is the effective number of neutrino species, which does not need to be an integer. With this parameterization, 
the total radiation energy density is 



P7 



7 



11 



4/3 



off' 



p^(l + 0.2271 ATeff). 



(14) 



While photons interact with baryons efficiently at z > 1090, neutrinos do not interact much at all for z <C 10^*^. As a 
result, one can treat neutrinos as free-streaming particles. Here, we also treat extra ra diation sp ecies as free-streaming. 



With this as sumption, one can use the measured Cj"^ spectrum to constrain A^eff |Hu et al. 1995; Hu et al . ^^ 
Bowen et al.||2002[ [Bashinsky fc Seljak||2004[ ). Section 6.2 of |Komatsu et al.| ('2009) ancTpTTot Komatsu^t^L (20T 
discus s pre vious attempts to c onstrain iVeg from the CMB and provide references. More recently, Dunklcy et aTI] 
(2011) and Keisler et al. (2011) constrain TVes using the seven-year WMAP data combined with ACT and SPT data, 
respectively. In this paper, we assume the sound spe e d and anisotropic s t ress of any extra radiation species are the 
same as for neutrinos. See Archidiacono et al. (2011); Smith et al. (2012); Archidiacono et al. (2012) for constraints 
on other cases. 

Neutrinos (and por) affect the power spectrum, Cj^ , in four ways 



Figure [8] compares models with A'cff 
possible. 



To illustrate and explain each of these effects, 
3.04 and iVcff ~ 7, adjusted in stages to match the two spectra as closely as 



Peak locations - The extra radiation density increases the early expansion rate via the Friedmann equation, 
_ &2^(^p^ _l_ result, increasing A'off from 3.04 to 7 reduces the comoving sound horizon, r^, at 

the decoupling epoch, from 146.8 Mpc to 130.2 Mpc. The expansion rate after matter-radiation equality is 
less affected, so the angular diameter distance to the decoupling epoch, dA, is only slightly reduced (by 2.5%). 
Therefore, increasing N^e reduces the angular size of the acoustic scale, 0* =rs/dA, which determines the peak 
positions. A change in 6^ can be absorbed by rescaling I by a constant factor. In the top-left panel of Figure [8j 
we have rescaled I for the A'efi = 7 model by a factor of 0.890, the ratio of 6** for these two models {9^, = 0?5961, 
0?5306 for Ncs — 3.04, 7, respectively). This rescaling brings the peak positions of these models into agreement, 
except for a small additive shift in peak positions; see [Bashinsky fc Seljak_ (2004] ). 

Early Integrated-Sachs- Wolfe effect - Extra radiation density delays the epoch of matter-r adiation equality 
and t hus enhances the first and second peaks via the Early Integrated-Sachs- Wolfe (ISW) effect ( Hu & Sugiyama 
1995). This effect can be compensated by increasing the cold dark matter density in the N^s = 7 model from 
TTji^ — 0.1107 to 0.1817, which brings the matter-radiation equality epoch back into agreement. (We do not 
change flbh'^, as that changes the first-to-second peak ratio.) The top-right panel of Figure [s] shows the spectra 
after making this adjustment. Note that changing ftch^ also changes 6**, so the / axis is rescaled by 0.957 for the 
Nes = 7 model in this panel. 

3. Anisotropic stress - Relativistic species that do not interact effectively with themselves or with other species 
cannot be described as a (perfect) fluid. As a result, the distribution function, /(x,p,t), of free-streaming 
particles has a non-negligible anisotropic stress. 



(2^ 



P PiP] 



/(x,p,t). 



(15) 



as well as higher-order moments. This term alters metric perturbations during the radiation era (via Einstein's 
field equations) and thus temperature fluctuations on scales I > 130, since those scales enter the horizon during 
the radiation era. On larger scales, fluctuations enter the horizon during the matter era and are less affected by 
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Fig. 8. — An illustration of four effects in the CMB anisotropy that can compensate for a change in the total radiation density, pr, 
parameterized here by an effec tive number of neutrino species, A^cfl- The filled circles with er rors show the nine-year WMAP data (in 
black), the ACT data (in green, [Pas et al.|2011[ |, and the SPT data (in violet, [Keisler et al.|20lT| |. The dashed lines show the best-fit model 
with N^ff = 3.04, while the solid lines snow models with N^ff = 7 with selected adjustments applied. (The other parameters in the dashed 
model are ^th^ = 0.02270, U^h'^ = 0.1107, Hq = 71.38 km/s/Mpc, = 0.969, A'^ = 2.384 X IQ-^, and r = 0.0856.) Top-left: the i-axis 
for the N^ff = 7 model has been scaled so that both models have the same angular diameter distance, d^, to the surface of last scattering. 
Top-right: the cold dark matter density, ^Ich^, has been adjusted in the N^g = 7 model so that both models have the same redshift of 
matter-radiation equality, Zcq- Bottom-left: the amplitude of the N^g = 7 model has been re- scaled to counteract the suppression of power 
that arises when the neutrino's anisotropic stress alters the metric perturbation. Bottom-right: the helium abundance, Yp, in the Ngg = 7 
model has been adjusted so that both models have the same diffusion damping scale. 



this term. Temperature fluctuations on these scales are given by the Sachs- Wolfe formula, ST/T = —TZ/b, while 



those on smalle r scales (ignoring the effect of baryons) are given by JT/T = - (1 + 4/^/15) '7^cos(fcr,) ( |Hu fc| 
Sugiyam"a||1996 1 , where fi, is the fraction of the radiation density that is free-streaming, 

MN..).^^^^^^. (16) 
j^y ctt; 1 + o.22717Voff 

The small-scale anisotropy is enhanced by a factor of 5(1 + 4fi,/15)^^ due to the decay of the gravitational 
potential at the horizon crossing during the radiation era. Since the anisotropic stress alters the gravitational 
potential (via the field equations), it also alters the degree to which the small-scale anisotropy is enhanced 
relative to the large-scale anisotropy. Therefore, the effect of anisotropic stress can be removed by multiplying 
CF^il ^ 130) by (1 + 4/^/15)2. In the bottom-left panel of Figure § we have multiplied Cf"^ at all I by 
[1 -h4/^(7)/15]V[l -h4/i,(3.04)/15]2, where U{7) = 0.6139 and /^(3.04) = 0.4084. The two models now agree 
well, but the iVofi = 7 model is greater than the standard model at I < 130 because the anisotropic stress term 
does not affect these multipoles. 

Enhanced damping tail - While the increased expansion rate reduces the sound horizon, rg, it also reduces 
the diffusion length, r^, that photons travel by random walk. The mean free path of a photon is Ac = l/(cTn-e)- 
Over the age of the universe, photons diffuse a distance w y^Sct/Xc \c oc \c/H , and fluctuations 
within rrf are exponentially suppressed (Silk dampin g, [Silkj|1968| ). Now, while the sound horizon is proportional 
to l/H, the diffusion length is proportional to due to the random walk nature of the diffusion, thus, 

Td/Tg oc \/H . As a result, increasing the expansion rate increases the diffu sion length relative to the sound 
horizon, which enhances the Silk damping of the small-scale anisotropy ( jBashinsky fc Seljak,,2004j ). Note that 
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TABLE 7 

Relativistic degrees of freedom and Big Bang Nucleosynthesis^ 



Parameter WMAP 


+eCMB 


+eCMB+BAO 


+cCMB+BAO+Ho 


Number of relativistic species'' 








iVeff > 1.7 (95% CL) 


3.89 ±0.67 


3.55 ±0.60 


3.84 ± 0.40 


Tis 0.988 ± 0.027 


q85+''-''l* 
u.9ao_Q Qj^g 


0.969 ±0.015 


0.975 ± 0.010 


Primordial helium abundance'' 








FHe < 0.42 (95% CL) 


0.299 ±0.027 


0.295 ±0.027 


0.299 ±0.027 


Us 0.973 ±0.016 


0.982 ±0.013 


0.973 ±0.011 


0.977 ±0.011 


Big bang nucleosynthesis'^ 








ris 


2.92 ±0.79 
302+0 038 

U.OUZ_(, (jgg 

0.978 ±0.019 


2.58 ±0.67 

r, nil +0.036 

u.oll_Q 

0.965 ±0.015 


n 978+0.034 

0.980 ±0.011 



^ A complete list of parameter values for these models may be found at 

|http: //la mbda . gsf c. nasa . gov/ 

" The parameters Ngff and Yhc comprise one additional parameter each in these table 
sections. 

The parameters N^f{ and Yhc a^r^ fit jointly in this section. 

rd/rs also depends on the mean free path of the photon, rd/vs oc yHTXc oc a/ H/ue, thus one can compensate 
for the increased expansion rate by i ncreasing the number density of free elect rons. One way to achieve this is 
to reduce the helium abundance, Yp (IBashi nsky fc Seljak|2004 Hou et ar]|2011 ): since helium recombines earlier 
than the epoch of photon decoupling, the number density ot t ree electrons at the decoupling epoch is given by 
(1 — Yp) Hf), where ni, is the number density of baryons (Hu et al. 



1995 



see also Section 4.8 of Komatsu 

et al. 2011). In the bottom-right panel of Figure [s] we show Cf ^ for the A^eff = 7 model after reducing Yp 
from 0.24 to 0.08308, which preserves the ratio rd/vs- The solid and dashed model curves now agree completely 
(except for I < 130 where our compensation for anisotropic stress was ad hoc). 

4.3.2. Measurements of N^g andYue: testing Big Bang nucleosynthesis 



Using the five-year WMAP data alone, Dunkley et al. (2009) measured the effect of anisotropic stress on the power 



spectrum and set a lower bound on A'cff . However, BAU and Hp data were sti ll required to set an upper bound due to 
a degener acy with the ma tter-radiation equality redshift ( Komatsu et al.|2009 ). This was unchanged for the seven-year 



analysis ([K omats u et al 
and AC 
nine-year 



20111. Now, with much improved me asurementsof the enhanced damping ta il from SPT 



1 ( ^2.2. ip, Civm "data alone are able to determine N^g (|Dunkley et al.||2011 

,r WMA 



Keisler et al. 2011). Using the 



WKP data combined with SPT and ACT, we find 

^eff = 3.89 ± 0.67 (68% CL) WMAP+eCMB; fixed. 

The inclusion of lensing in the eCMB likelihood helps this constraint because the primary CMB fiuctuations are still 
relatively insensitive to a combination of N^s and fim^^, as described above. CMB lensing data help constrain flmh^ 
by constraining erg. The measurement is further improved by including the BAO and Hq data, which reduces the 
degeneracy with the matter-radiation equality redshift. We find 

WMAP-HeCMB+BAO+iJo; >hc fixed. 

We thus find no evidence for the existence of extra 



iVeff = 3.84 ± 0.40 (68% CL) 



which is consistent with the standard model value of Ncg — 3.04 
radiation species. 

As noted above, this measurement of N^s relies on the damping tail measured by ACT and SPT, which is also 
affected by the primordial helium abundance, Ihc- Figure [9] shows the joint, marginalized constraints on Ncs and Yhc 
using the above two data combinations. As expected, these two parameters are anti-correlated when fit to CMB data 
alone (black contours). When BAO and Hq measurements are included, we find 



3.55 



+0.49 
-0.48 



F0.034 
-0.032 



(68% CL) 



WMAP+eCMB+BAO+iJn 



iVcff = 

Fhc = 0.278^ 

When combined with our measurement of the baryon de nsity, both of th ese values are within the 95% CL region of 
the standard Big Bang nucleosynthesis (BBN) prediction (l Steigman|2008 | , shown by the green curve in Figure |9] Our 
measurement provides strong support for the standard BBJN scenario, 'i'ablc [7] summarizes the nine- year measurements 
of iVoff and Fhg- 

4.4. Neutrino Mass 

The mean energy of a relativistic neutrino at the epoch of recombination is (E) = 0.58 eV. In order for the CMB 
power spectrum to be sensitive to a non-zero neutrino mass, at least one species of neutrino must have a mass in excess 
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Fig. 9. — Joint, marginalized constraints (68% and 95% CL) on the primordial helium abundance, Ihci and the energy density of 
"extra radiation species," parameterized as an effective number of neutrino species, A^eff . These constraints are derived from the nine-year 
WMAP+eCMB data (black), and from WMAP-|-eCMB-|-BAO-|-Ho data (red). The green curve shows the predicted dependence of Yhc 
on Ngg from Big Bang Nucleosynthesis; the dashed lines indicate the standard model: A^eff = 3.04, Yge = 0.248. 



of this mean energy. If one assumes that there are Ncs = 3.04 neutrino species with degenerate mass eigenstates, this 
would sug gest that the lowest to tal mass that could be detected with CMB data is '^rriv ~ 1.8 eV. Using a refined 
argument, Ichikawa et al. (2005) argue that one could reach ~1.5 eV. When we add 'Y^m^ as a parameter to the 
ACDM model we obtain the ht given in Table |8j specifically 

< 1.3 eV (95% CL) WMAP only, 

which is at the basic limit just presented. 

When the mass of individual neutrinos is less than 0.58 eV, the CMB power spectrum alone (excluding CMB lensing) 
cannot determine '^m^; however, tighter limits can be obtained by combining CMB data with BAO and Hq data. 
For a given ^ch^ and Hq, adding massive neutrinos results in a larger present-day total matter density, V,m, giving a 
smaller dark energy density for a flat universe, and hence a smaller angular diameter distance to the decoupling epoch. 
This change in distance can be compensated by lowering Hq. By the same token, for a given flch^+^i^h^ , adding 
massive neutrinos results in a smaller matter density at the decoupling epoch (since neutrinos are still relativistic 
then), which produces a larger sound horizon size at that epoch. Both of these effects cause the angular size of the 
acoustic scale, 9^,, to be larger, shifting the CMB peaks to larger angular scales. Furthermore, a reduced matter density 
at the decoupling epoch produces an earlier matter-radiation equality epoch giving a larger early ISW effect which, 
in turn, shifts the first peak position to a larger angular scale. This effect can again be compensated by lowering Hq. 
Therefore, indepen dent information on H q obtained from local distance indicators and from BAO data helps tighten 
the limit on ^ rrii, ( Ichikawa et al.|[2005 ). We find 



^ < 0.44 eV (95% CL) WMAP-f eCMB-hBAO+iJo, 



which is 25% lower than the bound of 0.58 eV that was set with the seven-year analysis ^Komatsu et al.|[2011 ). 

Since massive neutrinos have a large velocity dispersion, they cannot cluster on small scales. I'his means that a 
fraction of matter density in a low redshift universe (when neutrinos are non-relativistic) cannot cluster, which yields 
a shallower gravitational potential well, hence a lower value of as. As a result, one sees a clear negative correlation 
between us and J2 '"^f (s^e, e.g., the middle panel of Figure 17 of Komatsu et al.|20 09 |). Therefore, adding independent 
information on as obtained from, e.g^the abundance of galaxy clusters (Vikhlinin et al. 2009b Mantz et al. 20101 
helps tighten the limit on ^ m^. See §5] for a discussion of recent measurements of as from various cosmological probes 
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TABLE 8 
Neutrino mass^ 



Parameter 



WMAP 



heCMB 



hcCMB+BAO 



(-cCMB+BAO+Ho 



New parameter 

(eV) < 1.3 (95% CL) 
Related parameters 
as 0.706^ 

2.48 ±0.12 
0.962 ±0.016 



109 A^^ 



R+0.077 
"-0.076 
1 1 ^57+0-0048 

641+°'°^'5 



< 1.5 (95% CL) < 0.56 (95% CL) < 0.44 (95% CL) 



660+° °^® 
0.1183 ±0.0044 
0.586tO:0?o 
2.59 ±0.12 
0.947 ±0.014 



0.750 



,+0.044 
-0.042 



0.1133 ±0.0026 
0.695 ±0.013 

9 4^9+0.075 
^■^•^^-0.074 

0.9628 ±0.0086 



0.770 ±0.038 
0.1132 ±0.0025 
0.707 ± 0.011 
2.438 ±0.074 
9649+°-°°^^ 



A complete list of pa rameter values for this model may be found at 
[http :^/lambda ■ gsf c ■ nasa . gov/| 



TABLE 9 
Non-flat ACDM Constraints^ 



Parameter 



WMAP 



heCMB 



+eCMB+BAO +eCMB+_H'o +eCMB+BAO+J/() 



New parameter 














-0 037+°-0** 


-0.001 ± 0.012 




0.0049 ± 0.0047 


-0 0027+°-°°^9 


Related parameters 












Htot 


-1 „oy+0.042 
l.uo/ _Q Q44 


1.001 ± 0.012 


1 0049+°"°'"' 
^•'-"J*^-0.0041 


0.9951 ± 0.0047 


1 nn97+00038 




0.19 < n„ < 0.95 (95% CL) 


0.273 ± 0.049 


0.292 ±0.010 


0.252 ±0.017 


n 9oc:c:+0.0096 




0.22 < < 0.79 (95% CL) 


0.727 ± 0.038 


0.713 ±0.011 


0.743 ±0.015 


0.717 ± 0.011 


Ho (km/s/Mpc) 


38 < Ho < 84 (95% CL) 


71.2 ±6.5 


68.0 ± 1.0 




68.92i°:9^ 


to (Gyr) 


14.8 ± 1.5 


13.71 ± 0.65 


13.99 ±0.17 


13.46 ± 0.24 


13.88 ± 0.16 



^ A complete list of parameter values for this model may be found at http: //lambda .gsf c. nasa. gov/] 

such as cluster abundances, peculiar velocities, and gravitational lensing. 

4.5. Spatial Curvature 

The geometric degeneracy in the angular diameter distance to the surface of last scattering limits our abili ty to 



constrain spatial curvature, fifc, with primary CMB fluctuations alone ( Bond et al.|1997" Zaldarriaga et al]|1997 ). For 
example, the nine-year WMAP data gives a measurement with 4% uncertamty. 



-0.037 



-0.044 
-0.042 



WMAP-only, 



(see Table [9|). However, with the recent detection of CMB lensing in the high-^ power spectrum ( Das et aL]|2011 van 
Engelen et al. 20121, the degeneracy between ilm and JIa is now substantially reduced. This produces a significant 
detection of dark energy, and tight constraints on spatial curvature using only CMB data: when the SPT and ACT 
data, including the lensing constraints, are combined with nine-year WMAP data we find 



flA = 0.727 ±0.038 
ilk = -0.001 ±0.012 



WMAP+eCMB. 



Figure 43 in |Bennett et al. (2012) shows the joint constraints on {flm,^A) (and Qk, implicitly) from the currently 
available CMB data. Commning tne CMB data with lower-redshift distance indicators, such Hq, BAO, or supernovae 
further constrains Qk ( Spergel et al.|[2007 ). Assuming the dark energy is vacuum energy {w — —1), we find 



-0 0027+° °°3^ 



WMAP+eCMB+BAO+i/n 



which limits spatial curvature to be no more than .4% 

(12012). 



CL) of the critical density. These (flrn.^A) constraints 

are also shown in Figure 43 of [Bennett et al.[ 

The limits on curvature weaken slightly if dark energy is allowed to be dynamical, w ^ —1. However, with new 
distance measurements at somewhat higher redshift, where dynamical dark energy starts to become significant, the 
degradation factor is substantially less than it was in our previous analyses. We revisit this topic in §4.6[ 

4.6. Dark Energy 

The dark energy equation-of-state parameter, w = Pdo/pdc, where Pdc and p^c are the pressure and density of dark 
energy, respectively, governs whether pdc changes with time (w ^ —1) or not (w = —1). CMB data alone (excluding 
the effect of CMB gravitational lensing) are unable to determine w because dark energy only affects the CMB through 
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TABLE 10 

Dark energy constraints^ 



Parameter 



WMAP 



+cCMB +cCMB+BAO+/fo +eCMB+BAO+Ho+SNe 



Constant equation of state; flat universe 
w -1.71 <w < -0.34 (95% CL) 

Constant equation of state; non-flat universe 
w > -2.1 (95% CL) 

Qfe -0.052tO;O^i 

Non-constant equation of state; flat universe 
wo 



-1.07 



+0.38 
0.41 



-1.073! 



-1.19±0.12 

-0 0072+° ''°''2 

U.UU/Z_Q QQ43 



-1.34± 0.18 
0.85 ±0.47'' 



-1.037 



+0.071 
0.070 



-1.077 ± 0.072 
-0.0065 ±0.0040 



-1.17^ 



7+0.13 
-0.12 

M 
U.JO_o.4g 



^ A complete list of parameter values for these models may be found at Ihttp : //lambda . gsf c . nasa . gov/| 
^ The quoted error on Wa from WMAP+eCMB+BAO+i?o is smaller than that from 
WMAP+eCMB+BAO+Ho+SNe. This is due to the imposition of a hard prior, Wa < 0.2 - l.lioo, de- 
picted in Figure [To\ Without this prior, the upper limit on Wa for WMAP+eCMB+BAO+iTo would extend 
to larger values. 

1) the comoving angular diameter distance to the decouphng epoch, (i^(z*), and 2) the late-time ISW effect. The ISW 
effect has limited ability to constrain dark energy due to its large cosmic variance. The angular diameter distance to 
2* depends on several parameters (flm,^k, ^A, and Hq), thus a measurement of the angular diameter distance to 
a single redshift cannot distinguish these parameters. 

Distance measurements to multiple redshifts greatly improve the constraint on w. These include the Hubble constant, 
Hq, which determines the distance scale in the low-redshift universe; Z?y's from BAO measurements; and luminosity 
distances from high-redshift Type la supernovae. Gravitational lensing of the CMB also probes w by measuring the 
ratio of the angular diameter distance to the source plane (the decoupling epoch) and to the lens planes (matter 
fluctuations in the range oi z 1 — 2). Current CMB lensing data do not yet provide competitive constraints on w, 
though they do improve the CMB-only measur ement: co mpare the "WMAP^^ and "-|-eCMB" columns in Table [TO] 

New measurements of the BAO scale (i |2.2.2 ) and Hq (S 2.2.3) significantly tighten the 68% CL errors on a constant 
w for both flat and non-flat models 



_-| n7q+0 090 

i.U(0_Q Qgg 

-1.19±0.12 



(flat) 
(non-flat) 



WMAP-HeCMB+BAO+i?o- 



These constraints represent 35% and 55% improve ments, respectively, o ver those from the seven-year WMAP 
combination (see the fourth column in Table 4 of'Komatsu et a L]|2011 
(non-flat). Adding 472 Type la supernovae compiled by Conley et aL 



: w 



-1.10 ±0.14 (flat) and w 



BAO+Ha 
1.44 ±0.27 



(2011) improves these limits to 



-1.084 ±0.063 



-1.122 



+0.068 
-0.067 



(flat) 
(non-flat) 



WMAP+eCMB+BAO+i?o+SNe, 



where the errors inclu de systematic uncerta inties in the supernova data. Note that these limits are somewhat weaker 
than those reported in Komatsu et a l.| (2011 ), Table 4, c olumn 6, despite the smaller number of supernovae (397) in the 
"Constitution" sample compiled by Hicken et al. (2009), as that analysis did not include SNe systematic uncertainties 
in the seven-year analysis. 

Wa{l 



When w is allowed to vary with the scale factor according to w(a) 
Linder||2003[ ), we find, for a flat universq | 

+0.13 

WMAP+eCMB 

0.35+"-5° 



Wo 



a) (ChevalUer fc Polarski 2001 



Wo = -1.17] 



-BAO+Hn+SNe. 



-0.49 



10 shows the joint, marginalized constraint on wq and Wa- A cosmological constant {wo = — 1 and Wa = 0) is 
Doundary of the 68% CL region, indicating that the current data are consisten t with a time-independ ent dark 



Figure 
at the 

energy density. Comparing this measurement with the seven-year result in Figure 13 of Komatsu et al. (2011), we note 
that adding the new BAO and Hq data significantly reduces the allowed parameter space by eliminating < — 1. 

4.6.1. WMAP Nme-year Distance Posterior 

The "WMAP distance posterior" gives the likelihood of three variables: the acoustic scale, Ia, the shift parameter, 
R, and t he decoupling redshift, z^,. This likelihood is based on, and extends, t he original idea put forward by several 



authors (Wang & Mukherjee 2007 Wright 2007 Elgar0y & Multamaki 2007). It allows one to quickly evaluate the 



We consider only the flat case here since the non-flat case with Wa is not well-constrained by the present data. 
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-0.5 



W 







Fig. 10. — The joint, marginalized constraint on uiq and Wa, assuming a fiat universe. A cosmological constant {wq = —1, Wa = 0) is at 
the boundary of the 68% CL region allowed by theWMAP+eCMB+BAO+-ffo+SNe data, indicating that the current data are consistent 
with a non-evolving dark energy density. The shaded region is excluded by a hard prior, Wa < 0.2 — l.ltoo, in our fits. 



TABLE 11 

Inverse covariance matrix for 
the wmap distance posteriors 





Ia 


R 


2* 


Ia 


3.182 


18.253 


-1.429 


R 




11887.879 


-193.808 








4.556 



likelihood of various dark energy models given the WMAP data, without the need to run a full Markov Chain Monte 
Carlo exploration of the likelihood. 
Here, we provide an updated distance posterior based on the nine-year data. For details on how to use this simplified 



likelih ood, the definition of the above var i ables, and the limitation of this approach, see Section 5.5 of Komatsu et al. 



(|2011j) and Section 5.4 of K omatsu et al.| ( |2009[ ). 
'i'he likelihood is given by 

-21nL = (x-d)'^C~i(x-d), (17) 
where x = {I a, R, z.^,) are the parameter values for the proposed model, and the data vector d has components 

= l^MAP ^ 302.40 

d2 = i?^*'^^ - 1.7246 
is 



, ^ WMAP ^ ^ogo_j 



These are the maximum-likelihood values obtained from the nine-year data assuming a constant dark energy equation 
of state and non-zero spatial curvature (the 'OWCDM' model). The elements of the inverse covariance matrix, C~^, 
are given m Table [11] 

4.7. Constraints on Cosmological Birefringence 

If the polarization direction on the sky were uniformly rotated by an angle Aa, then some of the i?-mode polarization 
would be converted to i?-mode polarization. This can arise from a mis-calibration of the detector polarization angle. 
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but also from a phy sical mechanism called "cosmo logical birefringence," in which global parity symmetry is broken on 
cosmological scales ( Lue et ar]|1999 Carroll|[l998 ). Such an effect yields non- vanishing TB and EB correlations, hence 
non- vanishing Ur- A non-detection of these correlations limits Aa. 
A rotation of the polarization plane by an angle Aa gives the following transformation 



^TE,obs ^TE /OA \ 

C; =C; cos(2Aa), 


(18) 


^TB,obs r<TE ■ /OA ^ 
—(^i sm(2AQ;j, 


(19) 




(20) 


^BB,obs^^EEgi^2(2Aa), 


(21) 


Cf^^°'^^ = ^Q™sin(4Aa), 


(22) 



where the spectra on the right-hand side are the primordial power spectra in the absence of rotation, while the spectra 
on the left-hand side are what we would observe in the presence of rotation. We assume there is no _B-mode polarization 
in the absence of rotation. 

The low-Z TB and EB data at ^ < 23 yield Aa = -0?07 ± 4?82 (68% CL), while the high-Z TB data yield Aa = 
— 0?40 ± 1?30 (68% CL). The high-^ EB data are too noisy to yield a significant limit. Combining all the multipoles, 
we find 

Aa = -0?36± 1?24 (stat.) ± 1?5 (syst.) (68% CL). 
Here, we have added the s ystem atic uncertainty of ±1?5, to account for uncertainty in the WMAP detector polarization 



angle (Page et al. 2003 



statistical error has moc 



^2007 1. The WMAP limit on Aa is now dominated by this syste matic uncertainty. T he 
estly improved from 1?4 with the seven-year data (see Section 4.5 of Komatsu et al.||2011 ). 



5. OTHER CONSTRAINTS ON MATTER FLUCTUATIONS 



In this section, we summarize recent determinations of the matter fluctuation amplitude as traced by various measure- 
ments of large-scale structure. These include: cluster counts from optically-selected. X-ray-selected, and SZ-selected 
samples; measurements of N-point statistics in SZ maps, measurements of peculiar velocities, measurements of optical 
shear, and measurements of CMB lensing. To date, all of these observations are consistent with the WMAP nine-year 
ACDM fits, which give as = 0.821 ± 0.023 and as^^J = 0.434 ± 0.029. 

5.1. Cluster Observations 

Clusters are rare, high-mass peaks in the density field, he nce their number co unts provide an important probe of the 
matter fluctuation amplitude and, in turn, cosmology (see Allen et al. (2011) for a recent review). For cosmological 
studies, the main challenge with clusters is relating the astronomical observable (SZ decrement. X-ray flux, optical 
richness, etc.) to the mass of the cluster. Since the mass function is so steep, a small error in the zero-point of the 
mass-to-observable scaling can produce a significant error in the determination of cosmolo gical parameters. 

In the pas t two years, many new SZ-selected clusters have been reported by Planc k ( Planck Collaboration VIII 
2011| , ACT ( [Marriage et al.|2011[ penanteau et"al]|2012[ ), and SPT ( jWilhamson et al.J|20ll| jReichardt et al.||!j012b[ ), 
and they are providing new impetus for cosmological studies. Clusters are close to virial eqiiilibrium, so they should 
exhibit a tight relationship between integrated SZ d ecrement and mass; however, there are significant sources of non- 



Battaglia et al. 2012a,Cj ). Estimates of cluster 
i yZi measurements when one has both X-ray 



thermal pressure support that need to be modeled (Trac et al. 2011 

mass based on X-ray data agree well with estimates based on Planck SZ; measurements when one lias both A-ray 
and SZ data for the same cluster (Planck Collaboration X 2011 Planck Collaboration XI 20111, however, t here are 
intriguing discrepancies between some est imates based on optical and SZ data. I'his is seen in both the Planck (Planck 
Conaboratio n XII||2011[) and ACT data ( Hand et al.'2011; Schgal ct al. '20121, and several groups are exploring this 
discrepancy ( Anguio et al.||2012 Rozo et al. 2012; Biesiadziiiski et al. 2012) . 

The abundance of optically-selected clusters provided early, s trong evidence th at vacuum energy dominates the 
universe today. Based on the number density of ri c h clus ters. Fan et al. (19971 measured ag = 0.83 ± 0.15 and 
f2,„ = 0.3 ± 0.1. A recent analysis by Tinker et al. (20121, cornbimng two observables from the Sloan Digital Sky 
Survey: the galaxy two-point correlation function and the mass-t o-galaxy number ratio within clusters, found crgf^^^ — 
0.465±0.026, with ft,n = 0.29±0.03 and erg = 0.85±0.06. Zu et al. (2012) used weak lensing measurements to calibrate 



the masses of MaxBCG clusters in the SDSS data; they find asi^lm/0-S25)°-^°^ = 0.828 ± 0.049 

X-ray-selected cluster samples also provide constraints on the amplitude of matter fiuctuations. X-ray data allow 
one to b oth select the sample and calibrate its mass under the assumption of hydrostatic equilibrium. Vikhlinin et al.| 
(2009b| analyze cosmological parameter constraints from their Chandra cluster sample. Fitting ACDM parameters 
to the cluster counts, they find as{^M /0.25)'~'''^'' — 0.813 with a statistical error of ±0.012 and a systematic error of 



±0.02, due to absolute mass calibration uncertainty. 

With the new SZ-se lected cluster samples, groups are cali brating the SZ-decrement-to-mass scaling us ing wcak- 
lensing measurements ( Marrone et al . 2012; Miyatakc et al. 2012 High et al. 2012 Planck Collaboration III 2012); 
X-ray measurements (jBona mente et al.||2012( ^Benson et al.||20TT| ); and galaxy velocity dispersions ( Siton eFal.|i2012| ). 
Remarkably, these different groups are reporting fluctuation amplitudes that are consistent with the WMAP ACDM 
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fluctuation amplitude. For example, Benson et al. ( 2011) report i78(^m /0-25)° '^° — 0.785 ±0.037 based on an analysis 
of 18 SZ-selected clusters from the SP'l' survey, while |5eEgal et al.| ( |2011[ ) report = 0.821±0.044 and w = -1.05±0.20 
based on a joint analysis of WMAF seven-year data and 9 optically-confirmed SZ clusters. 

The n-point correlation function of SZ-selected clusters provides compleme ntary information to clus ter counts (the 
1-point function). The 2-point function is a potentially-powerful probe of erg (Komatsu & Scljak'2002); however, it is 
sensitive to the low-mass end of the Y { M') sca l ing relation for clusters, which is subject to astrophysical corrections 



(Shaw et al. 2010 Battaglia et al. 2010 2012a). Reichardt et al. (2012a I use simulations and observations to calibrate 



the SZ power spectrum (2-point function); they apply this to the SFl data to find ag = 0.807 ± 0.016. Higher-order 
correlation functions are less sensitive to low-mass cluste rs, so these moments are less affected by non-th ermal processes 
and more sensitive to the matter fluctuation ampl itude dHill fc Sherwin|2012 Bhattacharya et al.|2012 1 . Measurement 
of the 3-point function in the ACT SZ data ( ^ Wilson et al. 12012 ) yields cg = 0.78^Q g4. All of these measurements are 
consistent with the WMAP nine- year measurement of Ug = u.»21 ± 0.023, assuming ACDM. 



5.2. Peculiar Velocities 



ill (|2012|) combine redshi ft-space distortion data from BOSS (Reid et al. 
iggleZ ( Blake et al.||2011 1 , with local measurements of the peculiar veloci 



Hudson & 



Galaxy pecul iar velocities provide another independent probe of gravitational pot ential fluctuations. 

"20121, 6dFGS (IBeutler et al.||2012) 
iVTiild, to find ii™ = 0.i>59 i! 0M 



Turnbull 

and W 

(Tg = 0.748 ± 0'.U35, and a growi!h rate of 7 = rflnD/dlna = 0.619 ± 0.514. 
The kinematic Sunyaev-Zel'dovich (kSZ) effect can probe peculiar velocity fields over a wide range of redshift. 



Hand 



et al. (2012) report the first kSZ measurements of peculiar velocities at z ~ 0.35; they detect a signal consistent 
with predictions from N-body simulations that are based on the WMAP seven-year ACDM parameters. Future kSZ 
measurements should be able to provide precision tests of cosmology. 



5.3. Gravitational Lensing 

There are a number of complementary gravitational lensing techniques that measure the amplitude of potential 
fluctuations; 

CMB lensing - Large-scale structure along the line of sight deflects CMB photons and imparts a non-Gaussian 
pattern on the C MB fluctuation field. This is most easily detected on scales smaller than those probed by WMAP. 
Das et al. (2011 ) reported the first detection of CMB lensing using a measurement of the 4-point correlation function 
in the AC'i' ternperature maps; they report a deflection amplitude of A^ = 1.16 ± 0.29, w here A^ = 1 is the value 



predicted by the WMAP seven-year ACDM model. Using a similar technique on SPT data, van Engelen et al. (20121 



report = 0.90 ± 0.19. Since cx erg, this corresponds to an 8% measurement of Ug. When these CMB lensing 
measurements are combined with WAIAP s even-year data, they prov ide strong evidence for Dark Energy based purely 
on CMB observations ( .Sherwin et al. 2011] [van Engelen et al. 2012). 

Cos mic shear - Measu rements of cosmic shear in large optical siirveys directly probe matter fluctuations on small 
scales. Huff et al. (2011) analyzed 168 square degrees of co-added equatorial images from the Sloan Digital Sky 



Survey (SDSS) and found gg = .636±°;^°^ (when other cosmological parameters are fixed to the WMAP seven-year 



ACDM values). Lin et al. (|2012[) analyz ed 275 square degrees of co-added imaging from SDSS Stripe 82 and found 



o-g 



u.zoz,_g 052. 



Jee et al. 



(2012) report flm = 0.262 ± 0.051 and erg = 0.868 ± 0.071 from a cosmic shear 
study using the Deep Lensing Survey; wh en their results are co mbined with the WMAP seven-year data, they find 
J7m = 0.278±0.018 and erg — 0.815±0.020. Semboloni et al. (201l| analyzed both the second and third-order moments 
of the cosmic shear field in the HST COSMOS data; they found crg(f]A//0.3)° '*9 = 0.78^1^:26 using the 3-point statistic, 
in agreement with their result from the 2-point statistic: CTg(r27\f /0.3)°'^^ = 0.70^0 ■}4. 

Cross correlation - Correlating the cosmic shear field with the large-scale g alaxy distribution measures galaxy 
bias: the relationship been galaxies and dark matter. Mandelbaum et al. (2012) measured this cross-correlation in 
the SDSS DR7 data and used the inferred bias to determine cosmological param eters. They report erg( ^2^70.25)°-^^ — 



0.80 ± 0.05, where the errors include both statistical and systematic effects. Cacciato et al. (20121 use combined 



SDSS measurements of galaxy number counts, galaxy clustering, and galaxy-galaxy lensing, together with WMAP 
seven-year priors on the scalar spectral index, the Hubble parameter, and the baryon density, to find n„i = 0.278 

and as = 0.763t°:S5^^ (95% CL). 
Strong lensing - The statistics of lensed quasars prob e s the amplitude of fluctuations, the shape of galaxy halos 



4-0.023 
-0.026 



and the large-scale geometry of the universe. Oguri et al. (2012) have analyzed the final data from the SDSS Quasar 



Lens Search (19 lensed quasars selected from 50,836 candidates). They claim that the number of lensed quasar are 
consistent with predictions based on WMAP seven-year parameters. Assuming the velocity function of galaxies does 



not evolve with redshift, they report fl^ — 0.79_q 



+0.06 -1-0.06 
07 -0.06 



, where the errors are statistical and systematic, respectively. 



6. ON THE SZ EFFECT MEASURED BY WMAP AND Planck 



In Komatsu et al. 



([2011 1, we demonstrated that WMAP is capable of detecting and characterizing the Sunyaev- 
Zel'dovich (SZ) effect: tEe change i n CMB temperature due to inverse Compton scatterin g of CMB photons off 
hot electrons in clusters of galaxies (Zel'dovich & Sunyaev 1969 Sunyaev & Zel'dovich 1972 [). After o ur paper was 
published, the Planck collaboration published their first nieasurements of the SZ effect ( Planck Collaboration Vfll] 



24 




Fig. 11. — A compilation of the {flm, cs) constraints from large scale structure observations, discussed in qS] compared to the constraints 
obtained from CMB, BAO, and Hq data. The various large scale structure probes do not separately consErain the two parameters, and 
have somewhat different degeneracy slopes among t hem, but t hese independent measurements are quite consistent. The following la- 
regions are plotted: (a) crgQj^S = 0.465 ± 0. 026 from|Tinker et a l. (20121; (b) crs{n„i/0.325)°-^°^ = 0.828 ± 0.049 from|Zu et al. (201^; 
(c) o-8(f^m/0.25)'' *''' = 0.813 ± 0.032 from IVikhlinin et al.| ||2009bil ; (d) o-gCQm/0. 25)0-3 = 0.785 ± 0.037 fro m [Benson et al.| | |201l| ); 
(e) as(n,^/0.3)°" «~ -w.oo 1 11 1 



0.70 



+0.11 
0.14 



from 





pb^ 


(12012 






A^ilig I 



and 



Semboloni et al.l I 2011 1; (f) o-gf^m^ = 0.252 



from 



mmor axes are given by Sl„ 



^0.032 
-0.052 

= 0.259 ±0 



Lin et al, 
045 and 



(201211; (g) WMAP o nly; (h) 
0-8^1^748 ± 0.035 from [Hudson! 



o PJancJf's higher sensitivity and angular resolution, their measurements improve substantially upon 



In this section, we do not report any new results from the WMAP nine-year data, but we compare our seven- 
year findings with the correspondi ng Planck measurements. In addition, we no te that on- going blind SZ su rveys at 
arcminute angular scales by ACT (Marriage et al. 2011 Menanteau et al.||2012 | and SPT ( Song et ar]|2012 ) provide 



complementary information on clusters 

Coma cluster - Using the V- and W-band data, we are able to separate the SZ effect and the CMB fluctuation in 
the direction of the Coma cluster (Abell 1656). As a result, we find that the Coma cluster is sitting at the bottom 
of a —100 fiK CMB fluctuation, and that all the previous determinations of the SZ effect toward Coma that did 
not identify the primary CMB ov erestimated the SZ signal by about 25%. Our radial profile of the Coma cluster (see 
Figure 14 of|Komatsu e t al. | 2011[) is in ex cellent agreement with the much-improved radial profile measured by Planck 
(see Figure Abf Planck Collaboration X|[2012). 

The seven-year WMAP data are sensitive enough to measure 



Agreement -with X-ray-predicted SZ signal 

the SZ effect toward other nearby clusters. Among 49 z < 0.1 clusters with detailed Chandra observations (IVikhlinin 
et al. 2009a I, 29 are large enough to be resolved by WMAP and are outside the KQ75y7 sky mask. Among these, 
we detected the SZ effect in 20 clusters whose masses, M500, are greater than 2 x 10^"^ Mq. The Chandra data 
allow us to predict the SZ signal in each these clusters without relying o n any scaling relation s. We find very good 
agreement between the measured and predicted signals (see Figure 15 of Komatsu et al. 2011): when the Chandra- 
b ased prediction is fit to the SZ data from these 20 clusters, the best-fit amplitude is 0.82 ±0.12 (68% CL; see Table 12 



of Komatsu et al. 2011 1 . T his agreement has been confirm ed by the Planck collaboration with striking precision (see 
the left panel of Figure 4 of Planck Collaboration V|[2012 ). Their analysis used 62 2 < 0.5 clusters whose masses are 
greater than Afsoo = 2 x lO^'' Mq. 

Comparison -with the universal pressure profile - For a given cluster, the measured and predicted SZ signals 
agree well, if the prediction is derived from the detailed X-ray data on the same cluster. However, the agre ement is 
not as good if the prediction is derived from the so-called "universal pressure profile," (UPFp^ proposed by Arnaud 
et al.|(|2010|): the best-fit amplitude for the 20 clusters above M500 = 2 x 10^^ h'^ M© is 0.660 ± 0.095 (68% CL; see 



21 Here, the UPP refers to Equation (13) of 
= 0.12, and profile parameters C500 = 1.17' 



Arnaud et al 
, a = l.Ubl, /: 



(20101. Specifically, the pressure profile with a mass scaling of AjT^qq^"'' 
^^905, and 7 = 0.3081. 



/ith 
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Table 12 of Komatsu et al.pOlll. The Planck collabora tion observes the same trend with a similar magnitude (see 
the left panel of I'igure 4 of Planck Collaboration V||2012|). This may be caused by sample differences : the UPP is the 
median of a particular X-ray cluster sample, "KEXCJESH," derived from the ROSAT All-sky Survey ( Bohringer et al 



2007| . There is no guarantee that the median of the X-ray sample coincides with the median of our sample or the 
flanck sample. The disagreement between the UPP-based predictions and the WMAP SZ profiles means that WMAP 
is sensitive to details beyond average cluster properties, provided that the i nner structure of the cluster is re solved by 
the WMAP beam. The same is seen in the Planck analysis: see Figure 3 of Planck Collaboration XI (2011 1. 

Cool-core versus non-cool-core clusters - Motivated by this disagreement, we divided the samples into two 
sub-samples: (1) cooling-flow (or cool-core) clusters and (2) non-cooling-flow (or non-cool-core) clusters. Fitting the 
prediction from the UPP to the measured SZ data, we fin d best-fit amplitude s of 0.89 ± 0.15 and 0.48 ± 0.15 for sub- 
samples 1 and 2, respectively (68% CL; see Table 12 of Komatsu et al.pOll ). In other words, there is a statis tically 



significant difference between these two sub-samples. This is not so surprising: the X-ray data (Arnaud et al. 20101 



indicates that non-cooling-flow clusters have a significantly lower gas pressure in the core. We argued that this was 
the first time the same ef fect has been detected in the S Z data. The Planck collaboration has confirmed this (see the 
right panel of Figure 4 of [Planck Collaboration V|[20T2l) . 

Comparing to X-ray surveys - Combined, resolved measurements of the X-ray emission and SZ effect in a cluster 
give us a clear picture of the intra-cluster medium. For example, this has allowed us to detect the difference between 
cool-core and non-cool-core clusters in the SZ effect. However, such an investigation is limited to a small sample of 
clusters. Many more have been detected in the ROSAT All-sky Survey. How might we best use these clusters? 

In many cases, the only information available for these clusters is a redshift and an X-ray flux measured within a 
certain aperture. From this, one can derive an X-ray luminosity measured within a certain physical radius, for given 
cosmological parameters. We then need to use some scaling relations to relate the measured X-ray luminosity to the 
size (or the mass). This presents a challenge: while using more clusters increases statistics, using scaling relations 
introduces systematic errors. We tried three different scaling relations relating the X-ray luminosity, Lx, to the size. 



1 „ _ 0.753 h ^ Mpc r r /nrt-n 
1- ^500 - E(z) [^X/( W 



erg s 



0.228 



(Bohringer et al. 



2007 



derived from the Lx-temperature 



relation of Ikebe et al. ( [2002 ) and the size-temperature relation of Arnaud et al. ( 2005 ) , 



2. rsoo = 



3. rsoo = 



0.717 



£119(2) 



Mpc 



[Lx/{1Q^^ h-^ erg s-i)]°-222 ("REXCESS" scaling relation of 



Melin et al. 



2011), and 



0.745 h- 



£11^^(2) 



^[Lx/{1Q^^ h-^ erg s-i)]0-207 ("intrinsic" scaling relation of 



Mehn et al. 



2011) 



The mass (hence r^Qo) for the scaling relation 1 is estimated usin g the hydrostatic equ ilibrium, while that for the 



scaling relations 2 and 3 is estimated using the M^qo-Yx relation of Arnaud et al. (20071. Scaling relation 1 predicts 
the largest radius (hence mass) for a given luminosity, and the scaling relation 2 predicts the smallest radius for a given 
luminosity. Which relation should we use? There is no simple answer to this question, as different cluster catalogs 

have different selection functions^ 

For the study presented in |Komatsu et al. (20111, we use 499 clusters i n z < 0.2 whose Lx is greater than 
0.45 X 10"''* h^^ erg s~^. These clusters are taken from the REFLEX sample ( Bohringer et al.||2004 | in the southern 
hemisphere, and the extended BCS sample (Ebeling et al.|[l998 2000) in the northern hemisphere. We use the above 



scaling relations to relate Lx to 7*500, and convert it to the mass, M500. We then use this mass to calculate the 
predicted SZ signal from the UPP. By fitting the predicted SZ signals to the WMAP seven-year data, we find the 
best-fitting amplitudes of 0.59 ±0.07, 0.78 ±0.09, and 0.69 ±0.08 for the scaling relations 1, 2, and 3, respectively. For 
a sub-sample of high-Lx clusters with Lx > 4.5 x 10'*'' h-^ erg s"\ we find 0.67 ± 0.09, 0.90 ± 0.12, and 0.84 ± 0.11. 
Note that the quoted error bars are statistical-errors-on ly 68% CL. The unce rtainty in the scaling relation adds more 
errors to these numbers (see footnote (a) in Table 13 of Komatsu et al.|2 011| . 

This shows how uncertain this type of analysis is: depending on scaling relations one adopts, the results change 
significantly. For example, using the scaling relation 1, w e find that the UP P overpredicts the SZ effect. However, 
using the scaling relation 2 and the WMAP five-year data, Melin et al. (2011) find that the prediction from the UPP 
agrees with the measured SZ. Finally, the Planck collaboration uses the scaling relation 2 and finds that the prediction 
from the UPP agrees with the measured SZ. For the scaling relation 2, all of us agree (up to 2a): the predicted SZ 
from the UPP agrees with the measured SZ. However, for the other scaling relations the pr edicted SZ is rather low. 



The Planck collaboration finds the same result (see Figure 8 of Planck Collaboration X|2011 ). Therefore, these results 
indicate that the statistical analysis using many clusters and scaling relations is far less robust than that based on 
fewer clusters with detailed X-ray data. If the statistical analysis yields an unexpected result, then one should question 
the scaling relation. Thus, we do not feel there is any support for cosmological conclusions derived from stacking the 
WMAP data on the positions of hundreds of X-ray clusters, none of which is individually detected in SZ. Rather, one 
should conclude that a single simple X-ray based scaling law may not provide an adequately precise description of the 
wide variation in clusters. 

Note that the relatio ns 2 and 3 are derived originally for Lsoo, which came from the analysis of XMM-Newton observations of the 
ROSAT-det ected clusters ||Piffaretti et al.||2011| l, while we use the published values of Lx from the REFLEX and the extended BCS 
samples. In |Komatsu_e t aL T2PIT1 wc reported there was weak evidence indicating that lower mass clusters tended to be under-represented 
in their SZ signal. InjMclin et al. (2011J, this does not appear to be the case. We suspect the difference is due to the relation between Lx 
and L500 for low-lummosity clusters. 
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7. ACOUSTIC STRUCTURE IN THE 9- YEAR DATA 



7.1. Motivation 

In the standard model of cosmology based upon adiabatic scalar perturbations, temperature hot spots correspond 
to potential wells (i.e., over-dense regions) at the surface of last scattering; therefore, matter flows towards these hot 
spots. A crucial length scale in the CMB is the sound horizon at the epoch of decoupling, rs(z*); the angular size 
of the sound horizon sets the acoustic scale, 6 a = Tg/dA ~ 0?6. At twice the acoustic scale, the flow of matter is 
accelerating due to gravity, which creates a radial polarization pattern. At the acoustic sca le, the flow is decelera ting 
due to the central photon pressure, which creates a tangential pattern (see Section 2.2.2 of Komatsu et al.pOll and 
references therein). Around cold spots (potential hills), the polarization follows the opposite pattern, with tangential 
and radial polarization formed at 1?2 and 0?6, respectively. 

We detected this polarization pattern in the seven-year WMAP data with a statistical significance of Scr. The 



measured p attern was fully consistent with that predicted by the standard ACDM model (see Section 2.4 of Komatsu 
et al. 20111. Here we apply the same analysis to the nine-year data and find results that are again fully consistent 
with the standard model, now with a statistical significance of lOcr. 

The small-scale polarization data offer a powerful test of the standard model of cosmology. Once the cosmological 
parameters are determined by the temperature and large-scale polarization data, one can predict the polarization 
signal on small angular scales with no free parameters. This simple description is an important test of the standard 
cosmological model. 

7.2. Analysis Method 



Th e nine-year analysis replicates that of the seven-year data (see Section 2.3 and Appendix B of Komatsu et al. 
20111. We first smooth the foreground-reduced temperature maps from differencing assemblies VI through W4 to 



a common angular resolution of 0?5 (FWHM). We combine these maps with inverse- n oise- variance weighting, and 



remove the monopole from the region outside the KQ85y9 mask ( Bennett et al 
maxima and minima are obtained using the software hot spot in the MliiALFix dis 



20121). T he locations of the local 
irTBulion ( [Gorski et aT]|2005[ ). 

As in the seven-year analysis, we cull the hot spot list by removing all local peaks with T < (J, and vice versa. 
In the 66.34% of the sky outside the union of the KQ85y9 (temperature) mask and the P06 (polarization) mask, we 
find 11,536 hot spots and 11,752 cold spots remain in the temperature map. These counts are consistent with the 
expectation for a Gaussian random field drawn from the best-fit nine-year WMAP signal plus noise power spectrum. 

The raw polarization maps from differencing assemblies VI through W4 are combined using inverse-noise-variance 
weighting. We do not smooth the polarization maps for this analysis. We extract a 5° x 5° square region in the Stokes 
/, (5, and U maps centered on each hot and cold temperature spot. We combine the extracted temperature images with 
uniform weighting, while the Q and U images are combined with inverse-noise-variance weighting, excluding pixels 
masked by either analysis mask. Afterwards, we remove the monopole from the co-added Q and U images. There are 
625 0?2 pixels in each polarization image, which sets the number of degrees of freedom in the analysis. 

To make contact with the sta ndard model prediction, we work in a rotated polarization basis, Qr and Ur, introduced 
by Kamionkowski et al. (19971. These parameters are related to Q and U by 



Qr{0) = -Q(6/) cos 20 - J7(0) sin 20, 
Ur{0) ^ g(0)sin20- C/(6/)cos20, 



(23) 
(24) 



where 6 



'cost 



of Komatsu et al. (2011) 



is a position vector whose origin is the location of the temperature extremum; see Figure 1 
These Stokes parameters offer a simple test of the standard model, which predicts C/^ = 
everywhere, and Qr[t)} alternating between positive (radial polarization) and negative (tangential polarization) values. 

7.3. Results 

Figure [l2] shows the co-added T and Qr images from the nine-year data. We clearly see the alternating radial and 
tangential polarization pattern around the average hot spot, and vice- versa around the average cold spot. To test the 
agreement between data and theory, we fit the Qr maps to their predicted patterns, and let the amplitude be a free 
parameter. For the hot-spot Qr, we find a best- fit amplitude of 0.89 ± 0.14 (68% CL) with Ax^ = —41.3 relative 
to zero amplitude. For the cold-spot Qr, we find 1.06 ± 0.13 (68% CL) with Ax^ = —61.4. The combined best-fit 
amplitude is 0.973 ± 0.096 (68% CL). The data are fully consistent with the standard ACDM prediction, and the 
combined statistical significance of the detection is lOcr, compared to Sa for the seven-year data. 

The co-added Ur maps are consistent with zero. We fit the measured Ur maps to the predicted Qr patterns and find 
best-fit amplitudes of 0.02 ± 0.14 and 0.04 ± 0.13 (68% CL) around the average hot and cold spots, respectively. 



8. CONCLUSION 



We have used the final, nine- year W MAP temperatu re and polarization d ata (Bennett et al. 2012|) in conjunction 
with high-^ CMB power spectrum data ([Das et al.|2011 :; Keisler et al. 2011; R eichardt et al.|2U12a p3AO data (Bcutler 
et al. |20lTl Padmanabhan et a H2012 Anderson°eFarp012,^^Blake et ai.,2012|^ and a new Hq measurement (,Riess et aE 



to place strmgent constraints on the six parameters of the minimal ACDM model, and on parameters beyond the 
minimal set. The six-parameter model continues to describe all the data remarkably well, and we find no convincing 
evidence for deviations from this model: the geometry of the observable universe is fiat and dark energy is consistent 
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Fig. 12. — Co-added maps of temperature, T, and polarization, Q,., smoothed to a common resolution of 0?5, and stacked by the location 
of temperature extrema. (The polarization maps were not smoothed for the analysis, however.) Top-left: the average temperature hot 
spot. Top-right: the rotated polarization map, Qr, stacked around temperature hot spots. Bottom-left: the average temperature cold spot. 
Bottom-right: the rotated polarization map, Qr, stacked around temperature cold spots. The polarization images are color-coded so that 
the red {Qr > 0) shows the radial polarization pattern, while blue (Qr < 0) shows the tangential polarization pattern. The lines indicate 
polarization direction. These images are a striking illustration of BAO in the early plasma, and phase coherence in their initial conditions. 



with a cosmological constant. The amphtude of matter fluctuations derived from WMAP data alone, assuming the 
minimal model, erg = 0.821 ±0.023 (68% CL), is consistent with all the existing data on matter fluctuations, including 
cluster abundances, peculiar velocities, and gravitational lensing. The combined ( WMAP-|-eCMB-f BAO-f ffo) data 
set gives as = 0.820l^;[5i^ (68% CL). 

The basic predictions of single-field inflation models for properties of primordial curvature perturbations are well 
supported by the data: the te mperature fluctuations, which linearly trace primordial curvature perturbations, are 
Gaussian ( [Bennett et al. 2012) and adiabatic; they exhibit a slight power-law scale dependence, and the limits on 
primordial g ravitational waves are consistent with many inflation models, including one of the oldest, proposed by 
StarobinskyJ([W80) 



We find strong support for standard Big Bang nucleosynthesis from the joint constraint on the effective number of rel- 
ativistic species and the primordial helium abundance, which yields iVoff' — S.SSlp ^g and Ihc = 0-278^0 g32 (68% CL). 
The total mass of neutrinos is restricted to ^ < 0.44 eV (95% CL) without relying on information about the 
growth of structure 



We compared our seven-year measurements of the SZ effect, presented in [Komatsu et aL (2011), with the recent 
Planck measurements, finding that our results have been confirmed by Planck with striking precision. 

The improved polarization data around temperature extrema, now detected at lOcr, are in an excellent agreement 
with the prediction of the standard model based on adiabatic scalar fluctuations, providing a striking illustration 
of our physical understanding of the formation of acoustic waves in the early universe. No evidence for rotation of 
the polarization plane, e.g., by cosmological birefringence, is found: the nine-year WMAP bound is Aa = — 0?36 ± 
1?24 (stat.) ± 1?5 (syst.) (68% CL). The error is now dominated by systematic uncertainty. 

The nine-year WMAP data have reduced the allowable volume of the six-dimensional ACDM parameter space by 
a factor of 68,000 relative to pre- WMAP CMB measurements. When combined with the high-/ CMB, BAO, and Hq 
data the volume is reduced by an additional factor of 27. The maximum likelihood values of the ACDM parameters are 
given in Table [2] and the r nean and associated 6 8% CL error bars are given in Table [4] These results and those presented 
in the companion paper ( Bennett et al.|[2012 ) complete the WMAP Team's formal analysis and interpretation of the 
WMAP data. 
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APPENDIX 



PREDICTIONS OF STAROBINSKY'S i?^ INFLATION 
Primordial tilt 

As an example of a model that is consistent with th e WMAP nine- year data, we examine the predictions of Starobin- 
sky's inflation model. In 1980, |Starobinsky ( 1980 ) showed that one-loop quantum corrections to the Einstein-Hilbert 
action, which generate fourth-order derivative terms of 0{R^), where R is the Ricci scalar, lead to a de Sitter-type 
accelerated expansion of the universe. Starobinsky's motivatio n was not to solve the flatness and homogeneity prob- 



lems of the standard Big Bang model (this was later done by Guth [l981 ) , but to see whether one-loop corrections 
eliminate the classical singularity at the beginning of the universe. He was attempting to construct a cosmological 
model beginning from an initial de Sitter stage (not necessarily expanding) and ending in a radiation-dominated stage 
with a Friedmann-Robertson- Walker metric, with a mechanism for the graceful exit from inflation and the subsequent 
reheating phase. 

Starobinsky's work motivat ed Mukhanov and Chibisov in 1981 to consider quantum fluctuations in this model 
(Mukhanov & Chibisov 1981). They made the remarkable observation that quantum fluctuations generated during 
the de Sitter expansion are approximately scale invariant with a logarithmic dependence on wave number, and "could 
have lead to formation of gala xies and galactic clusters" (quoted from the abstract of their paper) . In current notation. 



Mukhanov & Chibisov 



(1981) show that A-jiik) oc (l -|- i In [see their Equation (9)], where H is the Hubble rate 



durmg milation and k is the comoving wave number. In the super-horizon limit, k <C aH, the primordial tilt observed 
in the wave number range accessible to WMAP, fcwMAP, is given by 



dlnA|(fc) 



d\nk 



fc — fcwMAP 



In 



fcwMAF 

aH 



2 



(Al) 



where N is the number of inflationary e-folds between the epoch when fcwMAP left the horizon and the end of inflation. 
With TV = 50, for example, one obtains = 0.96, whi ch is in agreement w ith our measurement. 

Among the one-loop correction terms considered by Starobinsky (1980), one particular term proportional to i?^ is 
found to be sufficient for driving inflation and creating curvature perturbations with the above spectrum, as the other 
terms vanish when the metric is conformally flat (e.g., a flat Friedmann-Roberts on- Wa lker metric). Then, one can 
obtain this result from the standard slow-roll calculation. As first shown by Whitt (19841, an action containing R and 
R^ is equivalent to an action containing R and a scalar field. Let us start with 



d'^xy/^ (R + aR^ 



(A2) 



Here, we have set 87rG — 1. By performing the conformal transformation, g^^i, — > = (1 -I- 2aR)g^^, and introducing 
a canonically normalized scalar field, ^E* = •\/3/2 ln(l + 2aR), Maeda^ (1988| ) finds that this system is described by a 
scalar field ^E* with a potential given by 



8a V 



1 



-y/2/3* 



(A3) 



which is quite flat for a large (also see Barrow fc: Cotsakis||1988 Salopek et al.||1989 ) 
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Once V{'^) is given, it is straightforward to compute the slow-roll parameters: 

2 yff 



(A4) 



where the primes denote derivatives with respect to ^. The number of e-folds of inflation between the end of inflation, 
te, and the epoch at which the observed perturbation left the horizon, t, is given, to the leading order in e, by 



N-: 



Hdt 



For the form of V{'^) given above, one finds, for a large value of ^Pe, 



-e 
3 



2/3* 



V 



4 

-e 
3 



2/3* 



N 



3 

-e 
4 



2/3* 



Therefore, the slow-roll parameters are given by N as 

3 



47V2' 



1 



i-e-, e < \v\- 

The primordial tilt is given by the slow-roll parameters as ( Liddle fc LythpOOO ) 



1 = -6e + 2r/ : 



2 



(A5) 



(A6) 



(A7) 



(A8) 



which agrees with the original result obtained by Mukhanov fc Chibisov (1981). 

Tensor-to-scalar ratio 



Prior to the invention of his inflationary model, Starobinsky ( 1979 1 calculated the energy spectrum of long- wavelength 
gravi tational waves produced during the de Sitter expansion, assuming the Einstein-Hilbert acti on (also see [G rishchuk| 
1975') . Such long-wave length gravitational waves induce tempe rature anisotropy in the CMB ( Rubakov et al.||1982f 
Fabbri fc Pollock 1983; iAbbott fc Wiselp84{ [Starobinsky|[l985| . 

Starobinsky (1983) calculated the energy spectrum of gravitational waves from his own inflation (also see ^ Mijic] 
[et al.,1986J . Instead of presenting their formulae, let us follow the standard slow- roll calculation. The tensor-to-scalar 
ratio, r, is given by the slow-roll parameter as r = 16e (Liddle fc Lyth 20001. Using the relation between e and N 
obtained in the previous subsection, one finds 

1 o 

(A9) 



12 



Therefore, while the inflation predicts a tilted power spectrum, the predicted tensor-to-scalar ratio is much smaller 
than 0(1 - n,) = 0{1/N) as it is of order 0{l/N^). 

Predictions of non-minimally-coupled inflation 
2/N and r — are exactly the same as those of an inflation model based on a 



These predictions, rig — 1 

scalar field nq n-minima lly coupled t o the Ricci scalar ( Spokoiny||1984 Accetta et al.|[l985 Futamase fc Maeda 1989 
Salopek et aL]ri989i |Fakir fc Unruhp90) 



J d^x^{l+C^^) R, 



(AlO) 



where the coupling constant, ^, is equal to —1/6 for the conformal coupling. For inflati on occurring in a large-field 
region, ^(j)^ 3> 1, with a scalar potential of V{(j)) cx (j>'^, the tilt is given by n.. - 1 = -2/N (|Salopek|l992||Kaiser|l995[ ), 
and the tensor-to-scalar ratio is given by (see Equation (5.1) of Komatsu fc Futamase||19g^ ) 

12 1 + 6£ , . , 

which becomes r 12/iV^ for ^ » 1 (also see Hwang fc Noh 1998"). 

This model has attracted renewed attention recently since ^Bezrukovfc Shaposhmkov| ( 2008 1 showed that the standard 
model Higgs field can drive inflation if the Higgs field is non-minimaliy coupled to the Ricci scalar with f ^ 1. The 
predictions of the original Higgs inflation for Us and r are therefore also Us 
eFaI][2009 ). 



1 = -2/7V and r = 12 /N^ (e.g., Bezrukov 
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